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pH ' Abstract 

This is the second of two papers in which we study real and virtual photon-proton 

scattering in a nonperturbative framework. In the first paper we have identi- 

p\ ' fied the leading contributions to this process at high energies and have derived 

expressions for them which take into account the renormalisation of the photon- 
quark-antiquark vertex. In the present paper we investigate the approximations 
and assumptions that are necessary to obtain the dipole model of high energy scat- 
tering from the results derived in the first paper. We discuss the gauge invariance 
of different contributions to the scattering amplitude and point out some subtleties 
related to gauge invariance in the correct definition of a perturbative photon wave 
function. As a phenomenological consequence of the dipole picture we derive a 
bound on the ratio of the cross sections for longitudinally and transversely polarised 
photons. This bound is independent of any particular model for the dipole-proton 
cross section and allows one to test the validity of the assumptions leading to the 
dipole picture in particular at low photon virtualities. We conclude that the naive 
dipole model formula should be supplemented by two additional terms which can 
potentially become large at small photon virtualities. 
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1 Introduction 

In this paper we continue our study of real and virtual photon-proton scattering in a 
nonperturbative framework which we have initiated in [1], hereafter referred to as I. 
Sections of that paper will be referred to as 1.2.1 and equations as (1.1) etc. We use the 
same notation and conventions as in I, and definitions made there will not be repeated 
here. In the present paper we include only references that are directly relevant to the 
material treated here, for a more extensive list of references we refer the reader to I. 

In the first paper we have classified the contributions to the process of real or virtual 
photon-proton scattering, see section 1.2.1, figure 1.2. We have identified the leading 
contributions at high energies in sect. 1.2.2 and have derived expressions for them 
which take into account the renormalisation of the photon-quark-antiquark vertex. 
In the present paper we study in detail the high energy limit. We investigate the 
approximations and assumptions that are necessary to derive the dipole model of high 
energy scattering from the results obtained in the first paper. We point out that the 
naive dipole model formula should be supplemented by two additional terms which can 
potentially become large at small photon virtualities. 

Gauge invariance is known to place strong constraints on current-induced scattering 
amplitudes. We find it therefore interesting to study the gauge invariance of different 
contributions to the scattering amplitude. It turns out that at finite energies the 
part of the Compton amplitude from which the dipole picture emerges is not gauge 
invariant separately. This also leads to subtleties in the definition of the wave function 
of longitudinally polarised photons which we discuss in some detail. 

As a key result of the present paper we identify the assumptions and approxima- 
tions that lead us from the most general description of the Compton amplitude to 
the usual dipole picture. We consider it a very important task to test the validity of 
those assumptions and approximations for the experimentally accessible ranges of the 
kinematic variables. In the dipole picture the total photon-proton cross section is a 
convolution of the square of the photon wave function with the reduced cross section 
describing the scattering of a colour dipole off the proton. At present, the reduced cross 
section cannot be derived from first principles and one has to retreat to models based 
on ideas like saturation. One finds that the currently available data permit considerable 
freedom for models of the reduced cross section. If one is to find phenomenological tests 
of the dipole picture (and hence of the assumptions on which it is based) one should 
thus concentrate on the other convolution factor in the dipole formula, that is on the 
photon wave function. This motivates us to study in some detail the properties of the 
wave function. 

As an immediate consequence of the dipole formula for the total virtual photon- 
proton cross section we derive a bound on the ratio R of the cross sections for longitu- 
dinally and transversely polarised photons that results only from the respective photon 
wave functions. Although there are only few data points available at energies suffi- 
ciently high for one reasonably to expect the dipole model to be applicable the bound 
turns out to give some indication as to the photon virtualities below which the dipole 
model becomes questionable. In this region the two additional terms mentioned above 
should become relevant. 

In the dipole picture one talks about a quark-antiquark pair scattering on a target. 
If we want to be rigorous we have to face the question how far we can go treating quarks 



as asymptotic states. Of course, quarks are confined and do not have a mass sliell. In 
this paper we shall nevertheless assume - as is usually done - that we can treat quarks 
and antiquarks as asymptotic particles. As we shall show we can then make a precise 
connection between the photon induced reaction and the reaction where the photon is 
replaced by a superposition of quark-antiquark states. We leave the investigation of the 
realistic case of quarks which have no mass shell for future work. At present we may 
give the following simple argument. Let us assume that the two-point functions for 
quarks have a Kallen-Lehmann representation with standard properties, but with no 
mass shell. This is indicated by calculations using lattice as well as Schwinger-Dyson 
equation methods, see for example [2]. We can then think that the mass-distribution 
function of a quark propagator may have a strong peak at some value mq which we shall 
call the quark mass. If the peak is approximated by a (5-function we have effectively a 
quark mass shell. 

This second paper is organised as follows. We start by discussing the high energy 
limit of the parts of the Compton amplitude that are leading at high energies and 
identify their asymptotic behaviour in section 2. In section 3 we investigate the separate 
gauge invariance of different parts of the amplitude. The derivation of the photon wave 
function is performed in section 4 where we also discuss the gauge non-invariance of this 
object and its consequences for the choice of the photon polarisation vectors. In section 
5 we generalise our findings to the case of more general amplitudes with incoming and 
outgoing photons and define dipole states. In section 6 we finally obtain the usual 
dipole picture and identify all underlying assumptions and approximations. Here we 
also study the properties of the photon wave function. In section 7 we derive a bound 
on the ratio of the cross sections for longitudinally and transversely polarised virtual 
photons that needs to be satisfied if the usual dipole picture is valid, and compare this 
bound to the available experimental data. Our conclusions are presented in section 
8. In appendix A we collect some technical details relevant to the discussion of the 
high energy limit in sections 2 and 4. Appendix B supplements section 3 with further 
discussion of the separate gauge invariance of different parts of the Compton amplitude 
which are subleading at high energies. 

2 The high energy limit 

In this section we shall study the high energy limit of the real or virtual Compton 
amplitude (1.4), |q| ^ oo with Q^ fixed. Our aim is to obtain eventually the usual 
dipole picture of high energy photon-proton scattering. As was discussed in section 
1.2.2 in the high energy limit the two amplitudes M^""' and M^^' give the leading 
contribution to the full amplitude. However, it is only the amplitude M^""' from which 
we expect the usual dipole picture to emerge. In the present section we will therefore 
study mainly that amplitude M^""' (1.14), (1.30). But it should be kept in mind that 
at low photon virtuality an important correction to the dipole picture can arise from 
the amplitude ^A^^' , see the discussion in section 1.2.2. 

We will mainly be concerned here with the term M^^'^' (1.43) of A^'") which we 
will find to be leading in the high energy limit, whereas the terms M^"'''^' to A^("'^) are 
subleading. In the amplitude A4'"'^^ by definition only the momentum k'^^' occurs but 
not the momenta k'^'^' to k'^^' which correspond to the other three parts of amplitude 



M^""' , see (I.B.26)-(I.B.29). Therefore we now simplify our notation and write k' instead 
of A;'(^) and k' instead of k'(i). 

Let us now consider the term 7W("'^) which we have obtained in (1.43) in the form 
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and let us in particular concentrate on the integration over lo in the complex w-plane. 
Recall that uj is the off-shell part of the energy of the quark on the left-hand side of the 
diagram in figure 1. 8a, see (LB. 7). We expect the matrix element of T^""' to be regular 
in the vicinity of w = 0. Thus we have the two explicit pole factors in the integrand 
(1) with the poles in a; situated at w = —ie and at a; = —/S.E + ie with 

A^ = fc° + k'^ - q^ , (2) 

see figure 1. One finds that always IS.E > 0. The poles will contribute a factor of order 
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Figure 1: The singularities of the integrand of M^^'^' (1) in the w-plane. There are 
poles at u; = —ie and uj = —AE + ie. 



1/q in the amplitude unless there is a pinching of the two singularities, that is for 

-A^ = g° - fc° - fc'° ^ . (3) 

The conditions for this to happen are easily derived and in principle well known. To 



make the article self-contained we discuss them here in detail. We have from (LB. 26) 
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where we suppose 



qkT = 0, (6) 

that is, longitudinal and transverse directions are defined relative to q. We then get 
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For large |q| the energy difference AE becomes small only if the terms proportional to 
q2 in the square brackets cancel. This happens for a 7^ 0, 1 if and only if 



that is for 
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Thus we get a pinch and a large contribution only if the splitting of the photon q is 
into a quark-antiquark pair where the longitudinal momenta are both in the direction 
q. Note again that we have used no arguments related to perturbation theory to arrive 
at this conclusion. 

We can rewrite (7) in yet another form: 
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Figure 2: The upper limit of k^ + 771,^ corresponding to m?{a,k.T) < Q^ (see (13), (14)). 



In figure 2 we show the range in k^ + rni corresponding to 



jS scale for 
< a < 1. We see that (12) corresponds to 



fl2) 



where Q is some chosen fixed mass scale for which we always suppose Q^ > Ami and 



kl + ml < 



Q' 



for a = and a = 1 , and to 



4(q2 + Q2) 
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for a = 1/2. Thus the pinch condition in the form 
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is only fulfilled for a ^ and 1 for very small k|, + mi whereas for q ~ 1/2 k^ + 
ranges up to 0{Q'^). This is discussed in more detail in appendix A.l. 

Taking now only the contribution from the pinch in (1) we obtain asymptotically 
for large |q| 
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Here the four-momenta k, k' are given as in (4) and (5) with < a < 1. We have 
inserted a ^-function to ensure the pinch condition in the form (15). This hmits k^ as 
shown in figure 2; see also the discussion in appendix A.l. 

We note that in A^("'^) to A^^"'^), (I.33)-(I.35), no pinch of the exphcit pole terms 
in Lo occurs and we expect these terms to be suppressed relative to the leading term 
(16) by a factor of the order 

AE^Q^ + Q^ 
qO - (^0)2 i ') 

for large |q|. 

In (16) we already see something like the dipole picture emerging. The photon 
splits into a qq pair with both the quark and the antiquark on the mass shell niq. This 
on-shell pair interacts with the proton to produce the final state. But the photon 'wave 
function' is not simply the vertex function. There is also the piece containing the kernel 
i^W:9 ) which in a way represents rescattering corrections of the qq-pair. 

We should keep in mind also the contribution from the amplitude Ai^^^ . In complete 
analogy to the discussion above we start from (1.45) and find a pinch singularity also 
for that amplitude. The leading term in M^ ' for large |q| then asymptotically becomes 
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The discussion of the high energy limit presented here is consistent with the well- 
known space-time picture of photon-hadron scattering at high energies that is the origin 
of the dipole picture. The quantity AE is the energy imbalance between the photon and 
the quark-antiquark pair, see (2). According to the energy-time uncertainty relation 
the quark-antiquark pair into which the photon fluctuates can therefore live at most for 
a time {AE)~^. In the high energy limit we have AE -^ (for fixed Q^) such that the 
maximal lifetime of the pair increases with increasing energy. In the rest frame of the 
proton, for example, the photon typically splits into the quark-antiquark pair already 
at a large distance from the proton. The actual interaction of the pair with the proton 
then takes place on much shorter timescales. Some typical values of AE in realistic 
situations are given in appendix A. 

In order to extract the dipole picture explicitly from the formulae obtained above 
it will turn out to be useful to address first the problem of gauge invariance of various 
parts of the amplitude. 

3 Gauge invariance 

In this section we want to study the question of electromagnetic gauge invariance of 
the real or virtual Compton amplitude. In particular we are interested in the question 
whether the different parts of the amplitude arising from the decompositions that we 



have performed up to this point are separately gauge invariant. In the following we 
will therefore consider the general amplitude without imposing the high energy limit. 
Only in section 3.3 we will come back to this limit. 

The full Compton amplitude (1.4) obviously exhibits electromagnetic gauge invari- 
ance. More precisely, electromagnetic current conservation implies 



quM'^r,{p',p,q) = o. 



(19) 



Here and in the following we always consider current conservation for the incoming 
photon, hence the contraction with q^,. The outgoing photon can be treated analogously. 



3.1 Skeleton decomposition of the full amplitude 

We first consider the general decomposition of the amplitude derived in section 1.2 and 
do not yet impose the high energy limit. 

Again we are in particular interested in the contributions M^""' and A^*-^' represented 
by figures 1. 2a and 1. 2b, which are not suppressed at high energies. Let us first consider 
the contribution M '-"^ , that is the term from which the usual dipole picture emerges at 
high energies, as we will see in the following sections. 

We start from expression (1.14) in which we recognise that the (/-dependence of 
Ai^"'' is fully contained in the term A^'^'^'^{q) given in (1.15). We therefore consider the 
expression 
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where we have integrated by parts. This contains the following expression, in which 
we can complete the derivatives to full covariant derivatives in order to obtain inverse 
quark propagators. 
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Inserting this back in (20) we obtain with (1.16) and (I. A. 5) 
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Consequently, the amplitude M.^"-' is separately gauge invariant, 

quM^^r{p',P,q) = ^- 



(22) 
(23) 



Next we turn to the amplitude M ^ ' . The factor in that amphtude which is relevant 
for our considerations here is according to (1.19) given by A^''^{q), see (I.C.3). Its 
contraction with q^, gives 
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where we have again integrated by parts. The derivative of the trace then becomes 
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Here we can again complete full covariant derivatives and obtain for (24) with (1.16) 
and (I.A.5) 
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Hence also the amplitude ^A^'^> is separately gauge invariant, 
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Finally we consider the remaining parts of the amplitude in the skeleton decompo- 
sition, A^(^) - Ai^^' , which are subleading in the high energy limit. Here we find that 
Ai^'^' taken together with Ai^ ' is gauge invariant, 



qu (-Mit'(p',P, q) + Mfrip'^P^ 9)) = , 



(28) 



whereas each of the remaining three amplitudes M^'^' , Ai^''' and M^^' is gauge invariant 
by itself, 

quM^;l''''{p',p,q) = qMPrip',P,q) = qMfr (p' ^P^Q) = «• (29) 

The derivation of these results is presented in appendix B. 

3.2 Decomposition of the amplitudes Ai^"-^ and tVI^'') 

Among the contributions to the skeleton decomposition of the full amplitude we have 
identified Ai^""' and Ai^^' as the ones that will be leading at high energies. We now 
turn to the spin sum decomposition which we have performed for these two amplitudes. 



At this point we still consider the general decomposition of these amplitudes performed 
before taking the high energy limit. 

We have decomposed ^A^"■> into the four terms ^A^'^'^\ j = 1, . . . , 4, see (I.30)-(I.35) 
and the graphical illustration in figure 1.8. According to (23) we clearly have 

q.t.^if'\p',P,q)=0. (30) 

i=i 

Since we have obtained the dipole picture only from M^"''^' , it is now interesting 
to see whether that part is gauge invariant by itself. We use the expression (1) for 
_y\^(a,i) g^j^(j remember that the terms in curly brackets in that equation equal Zq'y'^ via 

(1.41). Hence the contraction q^M^'g '^{p',P,q) contains as the relevant factor in the 
integrand the expression 

quUr{k)Zqj'^Vr'{k') = ZqUr{k) ^ Vr'^k') 

= ZqUr{k) [m-mq) + {^' + niq) - {^E)^"^] Vr'{k') 

= -{AE)ZgUr{k)-/^Vr'{k'), (31) 

where we have used (2), (4), and that by definition the spinors u and v are solutions 
of the free Dirac equation, 

Ur{k){}t — mq) = 0, 

{}t' + mq)Vr'{k') = 0. (32) 

Clearly, (31) is in general different from zero. Inserting that term back into the full 
expression (1) for M^""'^' we thus find that the resulting expression does in general not 
vanish. Hence the amplitude M^""'^' does in general not fulfil electromagnetic gauge 
invariance separately. 

Similar considerations hold for each of the other three terms M^""'^', j = 2, 3, 4, in 
(1.30). None of them is gauge invariant separately. 

The situation for the spin sum decomposition (I.C.ll) of the amplitude J\A^^' into 
the terms J\A^ '^' (j = 1, . . . , 4) is completely analogous. None of the terms J\A^ '^' given 
in (I.e. 13) - (I.e. 16) is gauge invariant by itself. In particular, the term M^ ''^' (I.C.14) 
gives rise to the same relevant factor (31) proportional to /S.E when contracted with 
the photon momentum. 

Finally, we turn to the decomposition of the photon vertex Zq^'^ in the amplitude 
_y\4(a,i) jj^^Q \^Q renormalised vertex term and the rescattering term, see (1.41) and the 
expression in curly brackets in (1). It is easy to see that none of the two contributions 
to M.^""'^' arising from these two terms is gauge invariant separately. Let us consider 
for example the renormalised vertex term. Contracting A^("'^) of (1) with q^ the term 
containing the renormalised vertex gives rise to a factor 

q^Ur{k)T^'i>{k,-k')Vr'{k') 

= Ur{k){k + k')^T^i>{k, -k')vr'{k') - (AS)n^(fe)r(«)°(fe, -k')Vr'{k') 

= ur{k) Wf~\k) - 4'^"\-fe')j Vr'{k') - {AE)Ur{k)T'^'''>\k, -k')Vr'{k') 

= -{AE)ur{k)T^'^^\k,-k')vr'ik'), (33) 



where we have used the QED Ward identity and the Dirac equation (32), and have 
assumed the existence of a mass sheh for quarks, that is Sp ("T-g) = 0. The remaining 
expression does not vanish and hence the corresponding term in M^""'^' is not gauge 
invariant separately. One can show that also the rescattering term is not gauge invariant 
separately. 

The same considerations hold for the analogous decomposition (1.41) of the ampli- 
tude >j(^'i) as applied in (1.45). 

3.3 High energy limit oi M^"-'^^ and M'^'"'^^ 

So far we have seen that the amplitudes Ai^"'' and M'^^' are gauge invariant by them- 
selves, but their most interesting parts in the high energy limit, Ai^^'^' and M^"'^' 
respectively, are not. Now we want to see what happens to the latter in the high 
energy limit. 

We start with the formula (16) for 7\/j("'i)^^ in the high energy limit. After contract- 
ing that amplitude with q^ the terms in curly brackets together with the two spinors 
give exactly the expression (31) proportional to AE, and this factor cancels against 
{AE)~^ in (16). Thus, q^Ai^'^'^'^^ has no (A£^)~-^ enhancement factor and is not of 
leading order at high energies. Consequently, we have for the ratio 

'■'^'°''""' o(^\, (34) 



|q|>l("'i)^o V |q| 

where we have chosen the component with z^ = in the denominator as a reference 
for the typical magnitude of the components of M^""'^'^^ . Since AE = 0(|q|~^) (see 
(A. 10)) this ratio tends to zero as |q| -^ co. This can be considered an approximate 
gauge invariance of M^'^'^'^'^ in the high energy limit in the sense that qyM^""'^'^^ is 
suppressed by two powers of |q| compared to the nominal power counting of the two 
factors of that contraction. 

It is a separate question whether q^M^"''^'^'^ itself vanishes in the high energy limit 
|q| -^ CO. That is in general not the case. We will see this explicitly when we discuss the 
proper choice for the polarisation vector in the photon wave function for longitudinally 
polarised photons in section 4.3 below. 

By the same argument as just given for A4*^"'^' also the gauge invariance of J\A^^'^' 
is restored asymptotically in the high energy limit, again in the sense that 



|q|_A4(6,l)/.0 \^ |q| y 

Finally, we turn to the decomposition (1.41) of the photon vertex Zq^^ into renor- 
malised vertex and rescattering term. In (33) above we have seen that the contraction 
of the renormalised vertex term with q^ gives rise to an expression proportional to 
AE. We can hence use the same argument as given above for the amplitude M^""'^': 
again, the {AE)~^ in (16) is cancelled by the factor AE from (33). It then immedi- 
ately follows from (1.41) and (31) that the same happens for the rescattering term. We 
conclude that approximate gauge invariance in the sense of (34) holds separately for 
the renormalised vertex term and for the rescattering term. The same conclusions hold 
also for the corresponding decomposition of the amplitude A^''''^'. 
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Note that the photon-quark-antiquark vertex Zq^'^ does not involve any momentum 
arguments. However, its separation (1.41) into renormahsed vertex and rescattering 
terms makes it necessary to assign momentum arguments for the quark Unes in the 
latter two. In (33) we have made use of the particular choice (1-42). Other choices 
for the momentum assignment would have given different results not necessarily giving 
rise to the separate approximate gauge invariance for the renormahsed vertex and the 
rescattering term in the high energy limit. Therefore the separate approximate gauge 
invariance of the two terms a posteriori justifies the momentum assignment (1.42). 

4 Photon wave function at leading order 

4.1 Definitions and general formulae 

In this section we stay in the high energy limit and derive the perturbative photon 
wave function at leading order from the results obtained above in (16) and (18). We 
will do this explicitly starting from the amplitude Ai^""' . Exactly the same steps can 
be performed also for M'^'''' . 

In the following we assume that the photon virtuality is sufficiently large so that 
we can apply perturbation theory. We work in the high energy limit and suppose 
furthermore 

Q |q| » Ik^l, mq , 
(1 — a) |q| » \kT\,'mq. (36) 

This explicitly excludes the end points in longitudinal momentum, a = and a = 1. 
In leading order of perturbation theory we have 

rii> = Y , (37) 

while the rescattering term vanishes in this order, 

i^^"?'"?') = . (38) 

We insert this in the leading contribution to the amplitude M^"'' at high energies as 
given in (16). Applying (5) we arrive at 

M, l^),p{p', s'W^^'^m', A', A'),q{k, A, A),p{p, s)) 
i;i%{a,kT,Q), (39) 



where we have defined 



i^%{a, kr, Q) = Qq ^^^^^^^po ^x(.k)Yvx'{k') 9 {Q^ - m\a, k^)) . (40) 

Note that we have separated the sums over r, r' into sums over Dirac indices A, A' and 
colour indices A, A' . In the following sections we will determine the leading contribution 
to this expression in the high energy limit. 
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We have introduced a normalisation factor N here which could in principle still be 
chosen freely at this stage. In order to conform with the usual normalisation of the 
photon wave function used in the literature we will pick the choice 

N = -2^/N~^ e y/ail^^ (41) 

with the proton charge e = ^/Atto^ and the number of colours Nc = 3. Note that the 
freedom in choosing A^ at this stage also involves a dependence on the fractions a and 
(1 — a) of the longitudinal momentum of the photon which are carried by the quark and 
antiquark, respectively. The particular dependence chosen here will be motivated later 
on by the way in which the remaining T-matrix element is related to dipole matrix 
elements and in which the latter depend on a and (1 — a), see sections 5.2 and 6.2 
below. 

Note further that (40) includes the theta function involving the mass scale Q intro- 
duced in section 2. We first notice that the function m?{a,\iT) becomes independent 
of q in the high energy limit |q| -^ oo, see (A. 12). Therefore this theta function natu- 
rally provides a regularisation of any singularities that might occur in the photon wave 
function at large transverse momenta |k7^|, or equivalently at small distances |Rt| in 
transverse position space. In the end of our calculation we want to choose the mo- 
mentum scale Q asymptotically large, thereby effectively removing the theta function. 
In view of this we will drop it already from now on in order to simplify the notation. 
However, it is very important to keep in mind that in our derivation of the photon wave 
function all possible singularities at small distances in position space are naturally reg- 
ularised. Therefore in our opinion the physical significance of any effects depending on 
the divergences of the photon wave function at small distances is doubtful. 

The Fourier transformation from transverse momentum to transverse coordinate 
space is given by 

i^\%,{a, Rt, Q) = I§^ ^^'"^" ^S'(«> kT, Q) . (42) 

We now define the photon wave function for transversely and longitudinally po- 
larised photons in momentum space as the contraction of V"! aa' ^^ (^^) "^i^h the re- 
spective photon polarisation vectors. For transversely polarised photons we define 

i,^^l%ia,kT,Q) = -e±u^%{a,kT,Q) , (43) 

and for longitudinally polarised photons we define 

i^^^ll, (a, kT,Q) = -eL ^Sa' ("' kr , Q) . (44) 

The corresponding wave functions ipl^^y{a,IlT,Q) and ipl^^y{a,IiT,Q) in coordinate 
space are obtained via the Fourier transformation (42). Suitable choices for the polar- 
isation vectors will be discussed in the following sections, where we give e±u explicitly 
in (51) and e^^^ in (57). Already here we should point out that the correct choice of 
a suitable polarisation vector in the case of longitudinally polarised photons involves 
some subtleties. As a consequence of that our definition (44) contains a polarisation 
vector e'j^ which differs from the most commonly used polarisation vector el, see (55) 
and (56) below. 
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For the discussion in the fohowing sections it is convenient to define the coordinate 
unit vectors in Minkowski space 



/ 1 \ 




/ \ 




/ \ 




/ \ 





, e(i) = 


1 


, e(2) = 





, e(3) = 





u 




U 




i 




u 


\o J 




\o J 




Vo^ 




V 1/ 



=(0) 



such that q points in the positive 3-direction, q = |q| e'^\ 
We choose conventions [3] for our spinors such that 



and 



u\{k) 



vxik) 



1 



y/FT 



nia 



(fc° + mg)xx 



o-kexA 



VfcO + mg V (^° + mq)exl 
where cr* are the Pauli matrices and 



1 
-1 



and for the 2-spinors xx we have the basis 



X+l 



X-1 

2 



It wiU be useful to introduce the quantity 



a(l — a)Q'^ + ml . 



(45) 



(46) 



(47) 



(48) 



(49) 



(50) 



4.2 Transversely polarised photons 

We start with transversely polarised photons, for which we choose the polarisation 
vectors 

/ \ 

V / 

Here the signs are chosen such that they are in agreement with the canonical Condon- 
Shortley sign conventions for angular momenta. 

With these polarisation vectors and with (40) the wave functions (43) for trans- 
versely polarised photons become 

i.%ia,kT,Q) = ^^Q^(^Ey \JkW ^A(fc)(7^ ±n^)^v(fcO • (52) 
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Using elementary algebra of Dirac and Pauli matrices and the expansions (A.7)-(A.ll) 
given in appendix A. 2 we obtain the leading contribution to this expression in the high 
energy limit as 



4^^^^%{a,kT,Q) 



1 



|±(/cri ± ikT2) a S^^^Sy-x - (1 - a) 5^^^i6x'-x 
+ mg6^^6x',x} ■ 



(53) 



Interestingly, terms proportional to |q| occur when the components of spinors and Dirac 
matrices are multiplied out in (52), but they cancel in the sum. We therefore have the 
important result that the photon wave function remains finite in the high energy limit. 
In order to obtain the wave function in transverse position space we apply the 
Fourier transformation (42) which gives with (A. 13) and (A. 14) 



^^'{t{a,RT,Q) 



^r. 



V2tt 






A,± 



.Sx' 



(l-a)6 



A,T 



6x' 



egKiiegRr) 



+ mq6)^^6x',xKoieqRT)j 



(54) 



where cpji = arg(i?n + iRT2) and Kj{x) are the modified Bessel functions. The results 
(53) and (54) are in agreement with those found in the literature up to the use of 
different phase conventions, see for example [4]. 



4.3 Longitudinally polarised photons 

Now we turn to longitudinally polarised photons. For these one often chooses the 
polarisation vector 



(el) 



i(|q|eW+gV=^) 



1 

Q 



( 1^1 \ 






(55) 



As we will explain momentarily, this choice is not appropriate for calculating the photon 
wave function in the high energy limit starting from the amplitude Ai^""'"^' found above. 
The correct result for the wave function of longitudinally polarised photons is obtained 
if one instead chooses the polarisation vector 









or explicitly in components 

(en ^ y^ 



,{0) 



,{3) 



Q 



|q| + 



Q' 



,{0) 



(56) 



.(3) 



1 

Q 



( H-q'\ 



\ 








(57) 



|q| / 
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With this choice for the polarisation vector and with (40) the wave function for the 
longitudinal photon (44) in transverse momentum space becomes 

^%{aM,Q) = -NQ,{AEr' 2Jk02k'0 ^^^ ^^^^^^^^ + 7')vx'{k') • (58) 
Using again the expansions (A.6)-(A.10) and 

q' = |q| + O (^) (59) 



|q| 

we obtain for the wave function for a longitudinally polarised photon in transverse 
momentum space 



r.(i)L (^, i_ ^^ _ o. /^ rr^n a{l - a)Q 

a(l -a)Q'^ + k^ + m2 



^;tA'(«'kT,Q) = -2Vn.V^Q, ,^ _^U^C2 ^^2 ^v,-a. (60) 



Hence the wave function remains finite for |q| ^ cx) also in the case of longitudinally 
polarised photons. In position space we obtain using (A. 13) 



ML 



^. 



ip: {y{a,'RT,Q) = V«em(3gQa(l " a) 6x' -X Ko{eqRT) , (61) 



vr 



again in agreement with the well-known result, see for example [4]. 

We now come to the problem of correctly choosing a polarisation vector for longitu- 
dinally polarised photons which we have already mentioned before. While the simplest 
possible choice for such a polarisation vector is e^ as given by (55) one can easily see 
that also 

^\ ^ |q| + /tg° (0) , g° + ^|q| (3) 



describes a longitudinally polarised photon for any real or complex number k. Indeed, 
since e'^ and e^ differ only by a multiple of q'^ one can invoke electromagnetic current 
conservation (19) to show that the contraction of the full Compton amplitude M^'^ 
with e'l^ gives the same result as the contraction with e^, 

e'LM':r,{p',p, q) = eLuM''J^{p',p, q) . (63) 

The same is true for the contraction of the two choices of a polarisation vector with 
the amplitude >J(")^'', 

sLM^l'^'ip^p, q) = eLMfr{p',P, q) , (64) 

because of (23). Due to (27) the same holds also for the amplitude 7\/j(^)^^. Therefore 
the polarisation vector e'^^^ is completely equivalent to elu as long as it is contracted 
with a gauge invariant amplitude, as expected on general grounds. The crucial point 
is now that in the definition of the wave function of longitudinally polarised photons 
(see (40) and (44) for the correct choice) the polarisation vector is contracted with an 
expression that is not gauge invariant, namely with u\{k)^'^v\i{k') which occurs in the 
amplitude JVl^'^'^>^^'^ . We have in fact seen explicitly in (31) that gauge invariance does 
not hold for this expression separately. 
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Let us see which concrete consequences this general observation has for the wave 
function of fongitudinally polarised photons. For that we calculate the analogue of (44) 
for eL. 

-e'LJ%i»M,Q), (65) 

which would naively appear to be a reasonable definition of the photon wave function 
for general k. With (40) and (62) we find for this contraction 

-NQ,{AE)-'^-^^ i uxik) [(|q| + Kg°) 7° " {q' + ^|q|) 7^] vx'ik') • (66) 

Now the product of the two spinors Ux{k) and v\i{k') with 7*^ occurring here gives the 
leading term 

^A(A:)7V(fc') = - } -7r= [2"(^ " "^'^ll' '^^''-^ + OKIq|, |kT||q|)] , 

(67) 
and interestingly their product with 7^ gives exactly the same leading term, 

nA(fc)7V'(^') = - "^ ^f7n== [2"(1 - «)|q|' 'Ja',-A + 0(m,|q|, |kT||q|)] , 

\J k + Vflq \J K + VClq 

(68) 
where we have again used (A.6)-(A.9). One finds in fact that even the terms propor- 
tional to ?Tig|q| and to Ik-pHql, that is the next-to-leading terms in the expansion in 
l/|q|, are identical in (67) and (68), but we will not need those terms here. For the 
square brackets and the spinors in (66) one then obtains 

(|q| + K g°) - (g° + K|q|) = (1 - /.) (|q| - (7°) (69) 

times the expression (67). Applying (A.6)-(A.10) and (59) and taking into account the 
prefactors in (66) we finally have 

««)"- <„, 1,^ ni _ _(, _ ,-!_ /^ rr^n a(y - a)Q 



-.J„«;V(«,k.,Q) = -(1 - „,yA-. v7S;:Q,-^_-i^^,_I|_,,,,_,^ (70) 

Obviously, this expression strongly depends on k which reflects the fact that W^^^y is 
not gauge invariant. Choosing arbitrary values for k here would lead to completely 
arbitrary wave functions for longitudinally polarised photons. 

What happens here is that cancellations take place between terms coming from 
different components in the contraction with the polarisation vector. If the chosen 
polarisation vector contains large terms proportional to |q| or g'', subleading terms in 
the amplitude (or in the product of spinors) can at the same time be enhanced due to 
these factors when the high energy limit |q| ^ 00 is taken. Those subleading terms, 
however, cannot be correctly determined only from V'i aa' '^^lich was extracted only 
from the amplitude A4("'^^ Instead, one would have to take into account also the 
subleading terms MS°-''^> to A^'"'^-*. Those terms do not have a pinch singularity (see 
the discussion in section 2) and are therefore difficult to identify in a nonperturbative 
framework using the methods that we have used here. We will comment on the situation 
in lowest order of perturbation theory at the end of the section. 
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If one wants to determine the correct wave function from the leading term at high 
energies alone, that is only from M^""'^' and hence from V'i aa" ^^^ must therefore 
choose a polarisation vector which remains finite in the high energy limit. Using (59) 
one sees that the only possible choice for k in (62) consistent with this requirement is 

n = -l + o(^), (71) 

up to subleading terms as indicated. The leading term in this expression turns (62) 
into the correct choice e'£ (56) all components of which remain finite in the high energy 
limit |q| — > oo. With that choice one obtains in fact the correct wave function (60) for 
longitudinally polarised photons as described above. 

Note that a similar problem with the choice of a polarisation vector would also 
occur for transversely polarised photons if one decided to add terms proportional to 
q'^ to e'^ (51). There, however, the simplest choices for the polarisation vectors do not 
involve any factors of |q| or q^. Therefore the gauge invariance problem is usually not 
discussed in the case of transversely polarised photons. 

Let us finally comment on the situation in lowest order of perturbation theory where 
the interaction of the quark-antiquark pair with the proton is given by two-gluon ex- 
change, possibly including the resummation of soft or collinear logarithms. Here the 
problem of choosing appropriate polarisation vectors can be studied explicitly, and the 
above results are found to be confirmed. In the perturbative situation it is possible to 
compute the photon wave function to lowest order starting from a gauge invariant ex- 
pression, as it has been done for example in [4]. There one starts with a set of diagrams 
in which the incoming and outgoing photons couple to a quark loop with two gluons 
attached in all possible ways. That amounts to four diagrams one of which is shown 
in figure 1.3. This set of diagrams is exactly the expansion of our amplitude J\A^"''^'^ 
(see figure 1. 2a) to lowest nontrivial order and constitutes a gauge invariant expression. 
For such a gauge invariant set of diagrams due to (64) and its perturbative expansion 
the different choices of photon polarisation vectors are fully equivalent. One can then 
compute the leading contribution to the quark loop at high energy and subsequently 
split the result into factors corresponding to the incoming and outgoing photon. More 
precisely, it is then possible to identify the perturbative photon wave functions from the 
different combinations of transversely and longitudinally polarised photons by invoking 
conservation of quark helicities at high energy. The resulting photon wave functions 
are in complete agreement with the results obtained above. 

5 Real and virtual photons in the initial and final states 

5.1 General formulae 

In this section we summarise our findings for the transition from photon to dipole 
scattering in the high energy limit. We give general formulae for real and virtual 
photons in the initial and final states. Let |a, in) be some initial, |6, out) some final 
state and let Ri{yi), . . . , RNivN) be some local operators. We consider the matrix 
element of figure 3, 

M''{q,yi,...,yN)= [ d^x e~''^^b,ont\T*Ri{yi)...RN{yN)r{x)\a,m). (72) 
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-f^*\q,iy) 




Figure 3: Diagram for the matrix element (72). 



Here T* is the covariantised T-product. In the following we will for the sake of brevity 
suppress the arguments yi and denote the above expression simply by Af^^q). For the 
case that no operators Ri are present M'^{q) is the amplitude for 



-f'^*\q) + a ^ b . 



As in previous sections we suppose 



-Q2<o. 



(73) 



(74) 



We should point out that in general we do not have qu-M-^ = for the amplitude 
(72) due to the T-ordering. Instead, q^M'^ vanishes only if for alH = 1, . . . , A'^ 



[j^ix),R,iy)]6{x^-y') = 0. 



(75) 



If the Ri are interpolating field operators for asymptotic states though, the full am- 
plitude will be described by the LSZ formula which then contains additional factors. 
If for example Ri corresponds to an outgoing particle with spin 1/2 and mass rrii, the 
LSZ formula will contain the additional integration 



(76) 



t / d^y,e^P^y^UsM)H^y.+m^)M-iq,...,y„...). 



The contraction of that full amplitude with qi, will then vanish. An example for this 
general mechanism is the separate gauge invariance of the amplitude M^-^' which is 
discussed in detail in appendix B. For the case that Ri is an electromagnetic current, 
on the other hand, the condition (75) holds in any case. Hence for interesting operators 
Ri current conservation for Ai'^ (or for the full scattering amplitude with LSZ factors 
for the case of interpolating fields Ri) is fulfilled. 

Our discussion in I and in section 2 can easily be adapted to the study of the high 
energy limit of the matrix element (72). We get for the leading term of Ai'^{q) for 
|q| ^ oo and fixed Q^ 



M%q) 



|q|^oo 






(fk 



1 



(27r)32fc02A:'0 A^ 
(6, out|T*i?i(yi) . . . RN{yN)\q{k' , r'),q{k, r), a, in) 
QqUr{k) [r(«)'^(A;, -k') + L'^''>{k, -k')] Vr'{k') . (77) 
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Here k, k' , AE are as in (2), (4), (5), and L^'^''^{k, —k') denotes the rescattering term as 
contained in (1.41), 

L(i>{k, -k') =Y,j K^''''''^S^pT^'>>S^p , (78) 

q' 
with the momentum assignments as in figure 1. 10. Note that the quark-antiquark state 
in the matrix element contains the appropriate renormahsation factor Z~^ as required 
by the LSZ formula as was made explicit in (I.36)-(I.39) for the case of the Compton 
amplitude. 

We have - of course having in mind the caveats pointed out in section 1 - also made 
the assumption that the quarks q have a mass shell. Furthermore a cutoff in transverse 
momenta is understood to ensure the pinch condition (15), see also the discussion 
concerning the corresponding theta function in section 4.1. 

Now we will use CPT = invariance to derive the analogue of (77) for real or virtual 
photons in the final state. Here and in the following we use the notation and general 
formulae for G-invariance as in [3]. With the antiunitary operator V{Q) representing 
a 0-transformation we have for the electromagnetic current 

{v{@)r{x)v-^{@))^ = -r{-x). (79) 

Let us define local operators Rf{y) by 

Rf{y) = {V{Q)Rj{-y)V-\e))^ (j = 1, . . . , iV) . (80) 

For state vectors we set 

y(G)|a,in) = I a®, out), 
y(e)|6,out) = |6®,in). (81) 

For - hypothetical - asymptotic quark and antiquark states we have 

l/(G)|g(/c,r),out) = -Y,\(i{k,t)M)Etr, 

t 

V{Q)\q{k,r), out) = ^\qik,t),m)Etr. (82) 

t 

Here r = (A,^), t = {t,B) denote again the combination of spin and colour indices as 
in I, and we have 

Etr = CrA ^BA ■ (83) 

Our starting point is (72), but with |a,in) replaced by |6 , in) and |6, out) by 
I a®, out). Furthermore we replace the Rjivj) by R^{y'A and take them in reversed 
order. This choice is possible since the states a and b and the operators Rj{yj) in the 
general formulae are arbitrary. We get for arbitrary points y'^, . . . ,y'j^ using (77) and 
(1.41) 

' d^x' e-^^"'(a®, out|T*i?® (y^) . . . i?f (y^) J^(x')|6®, in) 

d^k 1 



|q|-^oo 



^ E E / J^^^o^ Q,z,Mk)rvAk') (84) 



{a'',ont\T*R%{y'j,)...Rf{y[Mk',r'),q{k,r),b' 
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, m 



From (79) and (80) we have 

T*R%{y'j,) . . . Rf{y[) = {V{e)[T*Ri{y,) . . . RN{yN)]V~\Q))^ (85) 

and 

T*Rl,{y'i,)...Rf{y[).P{x') = {V{Q)[T*{-.P{x))Ri{yi) . . . RN{yN)]V-\@))^ , (86) 

where yj = — y' and x = —x'. Inserting this in (84) and using the rules of antiunitary 
operators we get 

d'^x e"^^{b,ont\T*{-J^'{x))Ri{yi) . . . RN{yN)\a,m) 

d^k 1 



> I 

q r',r t' ,t 



/d k 1 

(2vr)32fc02fc/0 A^ (-^t'r'Etr) Q,Z,Mkh'vAk') 

{q{k', t'), q{k, t), b, out|T*i?i(yi) . . . i?iv(t/iv)|a, in) . (87) 

Now we use 

'^EtrUr{k) = j5Vt{k), 
r 

J2Et'r'Vr'{k') = -j5Ut'{k') (88) 

r' 

for the Dirac times colour spinors, see (LB. 2) and (LB. 3), and get from (87) 
(i^xe*^"(6,out|T*(-J^(x))iii(yi)...i?jv(y^)|a,in) 

^ ^ ^ / i2j2k02k'0 XE Q'^^Mkh'uAk') 

q t',t '' ^ ' 

{q{k', t'), q{k, t), b, out\T*Ri{yi) . . . i?jv(y7v)|a, in) . (89) 



We interchange q and q in the final state, relabel the summation and integration vari- 
ables, and use (L41) in order to obtain 

d^x e*'?^(5,out|T*J^(x)i?i(yi) . . . i?7v(yjv)|a,in) 

d^k 1 



-^ I 



q r',r ^ ^ 

{q{k\ r'), q{k, r), 6, out|T*i2i(yi) . . . RN{VN)\a, in) . (90) 

This is our general result for the high energy limit of an amplitude with a real or virtual 
photon in the final state. 

5.2 Leading order formulae and dipole states 

In applications of the dipole formalism one frequently uses the photon wave function 
only to lowest order in a^. That is one makes the following replacements in (77) 

V^'i>{k,-k') -^ 7^ 

L^'i>{k,-k') -^ 0, (91) 
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and similarly in (90): 



L^^^^i-k^k) -^ 0. (92) 



This is justified for high enough Q^ where OsiQ'^) is small. 

Using now the definition of the lowest order photon wave function, see (40) and 
(41), we get from (77) 






<fk 



T 



q X',XA',A-"' " ^ ^ 



{b, out|T*i?i(yi) . . . RN{yN)\q{k , A , A ), q{k, A, A), a, in) 
1 1 

2^yl\cTT ^ya{l - a) ' 



6a'a -^==^=^= V^i'iv(«, kr, Q) . (93) 



With the Fourier transform (42) this reads 



eM^iq) > -i 



'll^"" q A',A 



EEi'^"/<*^«^ 



{b,ont\T*Ri{yi)...RN{yN)\D^''\q,a,-RT,X',X),a,m) 
^P%{a,RT,Q), (94) 

where we have defined dipole states 

|D('?)(q,a,RT,A',A)) = y f^e--'-^^^ ^ , ^ . 

' ^^' ' ^' ' ^^ Z^y (27r)2 2ViV;^^«(l-a) 

(5A'A|^(fe',A',A'),g(fc,A,A)). (95) 

Treating quark and antiquark as asymptotic states we find that these dipole states are 
eigenstates of three-momentum with eigenvalue q, 

P|D(^)(q,a,RT,A',A))=q|D(^)(q,a,RT,A',A)), (96) 

and satisfy for the kinematic region (36) the normalisation condition 

(D® (q, a, Rt, a', A)|d(«) (q, a, Rt, A', A)) = 6gg 5yy S-^^ilnf 2|q| S{a - a) 

(5(3)(q-q)(5(2)(RT-Rr). (97) 

A dipole state of the form (95) is for fixed R^ and a a superposition of quark-antiquark 
states with squared invariant masses (k + k')'^ for which we have at high energy due to 

(A.8) 

ik + kf >^T^,>^^,- (98) 

|q|-+oo a[i — a) a{l — a) 

Note that the mass of the quark-antiquark states is independent of Q^ in the high 
energy limit. The action of the squared four-momentum operator P^ on the dipole 
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states then gives in the high energy hmit |q| ^ co 
p2|D('?)(q,«,Rr,A',A)) = E / 

A I A 'J 



A' A' 



(Pkr it R 1 1 

J_ g-lkyRT ^^ - Sa> 

(27r)2 2ViVc^ V«(l - «) 



^^ Uq{k',X\A'),q{k,X,A)). (99) 

a(i — a) 

Note that also the mass of the dipole state is independent of Q^ at high energy. 

For photons in the final state we get a relation to dipole states analogous to (94) 
from (90) using the assumption (92): 



e I (fxe'^'' (6,out|T*J^(x)i?i(yi)...i?7v(yjv)|a,in) 



I'^l^"" , X',X 



^EE / ^" /^'^^ {'^i%i»,^T,Q)y (100) 

n A', A -^O -^ 

(Z)(«) (q, a, Kt, a', A), 6, out|T*i?i(yi) . . . i?^(y^)|a, in) . 



In the following section we will apply the general formulae (94) and (100) to the case 
of deep inelastic scattering. 

At this point a remark is in order concerning the interpretation of the dipole states 
introduced above. Frequently one thinks of dipole states as hadron states. However, 
there is still a small problem with the states (95) in this respect since they do not 
have the usual normalisation condition for hadrons. In particular, their normalisation 
depends on transverse position and the longitudinal momentum fraction of the quark 
and antiquark. But when considering the colour dipole as a hadron, these are continuous 
internal degrees of freedom which should not occur in the normalisation of the full 
hadronic state. This situation can be remedied by using dipole states smeared out in 
a and Ry. Consider smearing functions 

/r(Rr) > 0, 
fUa) > 0, (101) 

which are strongly peaked around Ry = and a = 0, respectively, and satisfy 

j (fRT\fT{^T)? = 1, 

(ia|/L(a)|2 = 1. (102) 

We can then define smeared dipole states as 

|D(«)(q,a,RT,A',A)) = I SRThi^T-^T) I dafL{a-a) (103) 

|D(^)(q,a,RT,A',A)). 
These states still satisfy (96), are normalised to 
(D® (q, a, Rt, a', A)|D(«)(q, a, Rt, A', A)) = 6^, 5^^ 6-^^{27rf 2|q| 5^^^ (q - q) , (104) 
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and can be thought of as hadron analogues. For the matrix element of the squared 
four-momentum operator between two such states we have 

(5®(q,«,RT,A',A)|p2|5(«)(q,a,RT,A',A)) 

= 6~yy 6-,^{27rf 2|q| 6^^H^ " q) {M^)a,^, ■ (105) 

Comparing with the normalisation (104) we see that {M'^)^-^ is the mean squared 
invariant mass of the smeared dipole state (103), and in the high energy limit we obtain 
using (99) 

(M')^,R. = I' da IfUa - a)p ^^^^ J ^ |/T(kT)P(k^ + m^) , (106) 

where we have defined the Fourier transform fx of the the smearing function Jt as 

/T(kr) = j <fRT e-^''^^^/T(RT) . (107) 

The latter satisfies with (102) 

/ ^ '^^^''^^'' = / ^'^^ 1/^(1^^)1' = 1 • (108) 

From (106) we conclude that at high energy the mean invariant mass of the smeared 
dipole states (103) depends only on a, niq and on the shape of the smearing functions 
Jt and f^, but is independent of Q^. 

6 Dipole picture for deep inelastic scattering 

6.1 Deep inelastic scattering 

As an important application of the general formulae (77) and (90) we consider the 
Compton amplitude as in (1.4) which we write now as follows 

{27:)H^*){p' + q'-p-q)M^r,{p\p,q) 

= -^— I (fx' (fx e^'i''''e-^'i''{p{p',s')\T*J^'{x')r{x)\p{p,s)) . (109) 
2vrmp J 

This amplitude should be understood as defined in the nonforward direction p' ^ p 
and the limit p' ^ p is understood in the following when we write M^'i'^{p^p^ q). That 
amounts to taking into account only the connected part of the matrix element in (109). 
The usual hadronic tensor of deep-inelastic electron-proton scattering is 

W'^''{p,q) = 1 J^l [A^7^(p,p,g) - {MZ{p,p,q)y] 6s>s (110) 

s' ,s 



nip \ q^ J \ q 
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with V = pq/nip. With Hand's convention [5] the transverse and longitudinal 7*-proton 
cross sections are [61 



aT{s,Q'^ 



2TTm> 



m. 



\e'i*e^W,,e\ 



2Tim> 



'^e^*e^W,..,.e'i 



mt 
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:iii) 
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27rm- 



V ^/^*„2 



27rm„ 
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(112) 



Here the 7*-polarisation vectors are as in (51), (55), (57). 

From (77) and (90) we get for the high energy limit of the Compton amplitude 

(2l)'M''l(p' + «' - p - «)>!"■' (p',p, q) 



TL 



S>1 



rp^^EE 



|q|,|q'Hoo 27rm„ ^^J (27r)32/02/'0 AE' ^^] (2TT)^2k02k'° AE 

' ' ' ' ' g' t',t -^ ^ ' q r',r "^ ^ ^ 



d^k 



Qq'Vt'il') [r(«')^(-/',/) + L(«')^(-/',/)] ut{l) 

{q'{l', t'),q{l, t),p{p', s'),out\ q{k', r'),q{k, r),p{p, s), in) 

QgUr{k) Ir^^^'^ik, -k') + L(«)^(fe, -A;')] Vr'{k') . (113) 

Here Ai?, k and k' are as in (2), (4), and similarly we define 



/ 



F + m^, 



I' 



q^ - 1)2 + m^, 
q'-l 



f + /'O 



,'0 



AE' = r + L"-q- 

Treating quarks and antiquarks in (113) as asymptotic states we can write 

{q'{l', t'), q'{l, t),p{p', s'),out\ q{k', r'),q{k, r),p{p, s), in) 
= 6fi + i{2TT)*6^'^\l' + l+p' -k' -k-p) 

{q{l', t'),q{l, t),p{p, s')\T\q{k', r'), q{k, r),p{p, s)) , 



(114) 



(115) 



where T represents the full T-matrix element for this scattering process with three 
incoming and three outgoing particles. The contribution 6fi gives rise to the discon- 
nected part of the matrix element in (109). It does not contribute to the cross section 
and will not be considered further here. In the high energy limit we have 



6^^Hl' + l+p' -k-k' -p) > 6^*\p' + q' - p - q) 



(116) 
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neglecting the terms AE and AE' in the energy sum on the l.h.s. With this we get 
from (113) 

^s'siP^P^l) |q|Jq,|^oo' 2^^^ J [271^21^21'^ AE>^^ J (2TT)^2k'^2k'0 AE 

{q'{l',t'),q\l,t),p{p',s')\r\q{k'y),q{k,r),p{p,s))\^^^^^ 
QqUrik) \T^'i>{k, -k') + L^'i>{k, -k')] Vr'{k') . (117) 

Note that when taking the high energy hmit of the Compton amplitude we want to 
keep only the contributions to the T-matrix element of (115) which are leading at high 
energies, as indicated by the subscript. In section 1.2.2 we have considered the high 
energy limit |q| ^ oo for the incoming photon and have found that the leading contri- 
butions are contained in the amplitudes Ai'^"'' and Ai'^^' of our skeleton decomposition 
of the Compton amplitude. In both of these amplitudes the incoming and the outgoing 
photon can be treated in a completely symmetric way. That symmetry can be seen 
from the diagrams representing these classes in figure 1.2 and has also been indicated 
at the end of appendix LA for the example of M^"'' . It is therefore straightforward 
to establish that the leading contribution to the matrix element in (117) is contained 
in the diagram classes corresponding to the amplitudes M^^> and A4' •* also for the 
simultaneous limit |q|, |q'| -^ oo. 

The result (117) is valid for forward as well as non- forward and real as well as virtual 
Compton scattering. 

6.2 The usual dipole picture 

To finally arrive at the standard formulae for the dipole model we consider forward 
Compton scattering. We need to make several assumptions (i)-(v) which we now discuss 
in detail. 

We start with two assumptions which we have already used before: 

(i) Quarks of flavour q have a mass shell niq and can be considered as asymptotic 
states. 

(ii) The rescattering terms L*-"^ '^ and L^'^>'^ in (117) are dropped and the vertex 
functions F^'^ •'^ and T^'^''^ are replaced by the lowest order terms in perturbation 
theory, that is ?(«')'' -^ jf^^T^i> -^ Y- 

We have made use of assumption (i) when we defined dipole states in section 5 above. 
Assumption (ii) has been used for deriving the photon wave function at leading order 
in section 4, and was applied also when we obtained the corresponding formulae for 
photons in the final state in section 5.2. Since the rescattering terms L^'^> start at order 
Ug in perturbation theory dropping them is consistent with keeping only the leading 
order contribution to the vertex function F''^'. In higher orders there will be correction 
terms to both the vertex functions and the rescattering terms. 
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With these assumptions we find from (117) with the photon wave function to lowest 
order as given in (40), (42) and with (94) and (100) 



e'M>:rAp,P^Q) 



asympt 



1 



27rm. 



P - - - JO 
^ Q A', A 



^Y. da SRrY.T. 



Q A', A 



da I d"i2T(^®^,(a,RT,Q) 



(Z)®(q,a,RT,A',A),p(p,s')|r|i?('?)(q,a,Rr,A',A),p(p,s)) 

7,AA' 



^P%ia,KT,Q), 



(118) 



where the subscript on the l.h.s. stands for the high energy limit, |q| -^ oo, and D 
are the dipole states defined in (95). At this stage only the leading contributions in 
the high energy limit should be included in the T-matrix element. As pointed out in 
section 6.1 above they are contained in the amplitudes A4'"' and Ai^ ' corresponding 
to the classes (a) and (b) of our skeleton decomposition of the Compton amplitude 
introduced in I. These two classes are shown again in figure 4, now for incoming and 
outgoing dipole states. The shaded blobs indicate again the functional integration over 



i:»®(q,a,RT,A',A) 




D(9)(q,a,RT,A',A) 



(b) 



i?(9)(q,a,RT,A',A) 




2?(«)(q,a,RT,A',A) 



Figure 4: The diagram classes containing the leading contribution to the T-matrix 
element of dipole states in (118) in the high energy limit. The classes (a) and (b) result 
from the amplitudes Ai^^' and Ai^ ' of figure 1.2, respectively. 



all gluon potentials with a functional measure including the fermion determinant. The 
discussion of typical diagrams and of the relative size of the diagram classes presented 
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in section 1.2.2 for incoming and outgoing photons holds in exactly the same way here 
for dipole states. In the following we will generically denote the two contributions to 
the T-matrix for dipole states corresponding to the two diagram classes of figure 4 by 
T^"-' and T'^', respectively. 

The following assumptions will concern the T-matrix element of dipole states. Since 
the dipole picture is usually applied only to the scattering of photons off unpolarised 
protons we will in fact need to make these assumptions only for the proton spin averaged 
T-matrix element 

^ Y, Ss'siD^'Hq, a, Rr, A', ~X),p{p, s')\T\D^'^\ci, a, Rt, A', X),p{p, s)) . (119) 

s' ,s 

We will nevertheless formulate the next two assumptions for the T-matrix element 
before averaging and point out some interesting issues concerning the proton spin in 
our discussion. 

Our third assumption is the following: 

(iii) The T-matrix element occurring in (118) is diagonal in the quark flavour, in a 
and Rt, and is proportional to the unit matrix in the space of spin orientations 
of the quark and antiquark in the dipole. 

As we will see this assumption is crucial for the dipole picture, but it is far from obvious 
whether it can be derived on general grounds. 

Our assumption (iii) is known to hold for the simplest perturbative diagrams (one 
of them is shown in figure 1.3), both for the case without or with the resummation of 
leading logarithms of the energy, see for example [7] . These diagrams also illustrate well 
the intuitive picture of high energy scattering on which the more general assumption 
(iii) is based. According to this picture a highly energetic quark or antiquark follows 
a straight trajectory and does not change its transverse position or its longitudinal 
momentum by a sizable amount while being scattered in the very forward direction, 
and similarly one expects its helicity to be conserved. For a discussion of this in a non- 
perturbative framework see for instance [8, 9]. The interaction would correspondingly 
be diagonal in flavour, in a and Ry, and in the helicity. This picture has been quite 
successful phenomenologically; it for example forms the basis of all eikonal formulations 
of high energy scattering. A fully nonperturbative derivation of this picture including 
an estimate of possible corrections could presumably be done along the lines of [8, 9], 
but this remains to be worked out. 

Let us for instance consider the diagonality of the interaction in the transverse po- 
sition of the quark. The idea that hadronic states in which the partons have flxed 
transverse positions become eigenstates of the T-matrix in the high energy limit has 
been put forward already in [10]. In the particular case of the dipole picture one of- 
ten argues that the typical lifetime of the quark-antiquark pair into which the photon 
fluctuates is much longer than the typical timescale of the interaction with the proton 
and that therefore the transverse positions (or other properties) of the quark and an- 
tiquark should not change during that interaction. But this argument does not make 
any reference to the actual dynamics of the interaction and in our opinion would need 
more rigorous support. While for soft gluon exchanges that argument can be made 
more concrete [8] the situation remains to be clarifled for general interactions. Interest- 
ingly, the results found in [11] indicate that the diagonality in Ry is violated already 
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in perturbation theory if one goes beyond the leading logarithmic approximation by 
taking into account exact gluon kinematics in the framework of /cT^-factorisation [12]. 
It would be very interesting to study this issue in more detail in the framework of a full 
calculation in next-to-leading logarithmic approximation. This could lead to a more 
solid test of the validity of assumption (iii) and hence of the foundation of the dipole 
picture at least in perturbation theory. 

An interesting issue is also the assumption that the helicity of the quark and the 
antiquark are conserved in the interaction. This helicity conservation has already been 
found in the perturbative study of various QED processes at high energy, see for ex- 
ample [13, 14, 15]. The results obtained there for the leading terms at high energy 
can be applied also for gluon exchanges between quarks and hence also hold in QCD. 
In perturbation theory one finds that there are also contributions involving a helicity 
flip if one goes to the next-to-leading logarithmic approximation, see for example [16] 
and references therein. It hence remains an open question to which accuracy the con- 
servation of quark helicities holds in general. An interesting possibility to address this 
question might be the study of spin-dependent structure functions as we will explain 
in section 6.3. 

As we have pointed out the intuitive picture of high energy scattering described 
above is reflected in the simplest perturbative diagrams, see for example figure 1.3. But 
one should remember that those diagrams contribute only to the amplitude M^""' and 
hence to the part T^""' of the full T-matrix element, see figure 4a. One can easily see 
that the diagonality in assumption (iii) is plausible only for that part T'"\ Let us 
consider the second contribution T^^' to the T-matrix element (corresponding to the 
part A^(^) of the amplitude), see figure 4b. One of the simplest diagrams contributing 
to T^ > is shown in figure 1.4. The flavours in the two quark loops in this diagram 
are independent of each other, and the contribution of identical quark flavours in the 
loops constitutes only the smaller part of the amplitude. Further it appears extremely 
unlikely that the outgoing quark and antiquark should have a high probability to be 
produced at the same positions in transverse space and to carry the same longitudinal 
momentum fractions as the incoming quark and antiquark. Similar considerations can 
be applied to the full matrix element represented by figure 4b. Hence we can at most 
expect that the part of the T-matrix element corresponding to the amplitude Ai''"'' 
satisfies the assumption (iii). 

In the light of this discussion of the simplest perturbative diagrams corresponding to 
T*-*' we can say that the above assumption (iii) is likely to exclude most contributions to 
the scattering amplitude contained in the amplitude A4' •*. However, we cannot make 
this statement rigorous based on fully nonperturbative considerations. Thus, although 
it is most probably a consequence of assumption (iii) we formulate it as a separate 
assumption of the dipole model: 

(iv) In the T-matrix element in (118) only the contribution T^"-* is kept while T' ' is 
neglected. 

Recall that the amplitude M^^' contributes only at higher orders of a^ compared to 
M-^"'' , but at low photon virtualities its contribution can well be important. 
Finally we need to make the following assumption. 

(v) The proton spin averaged reduced matrix element depends only on R|, and on 
s = {p + g)^. 
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This assumption implies in particular that the reduced matrix element is independent 
of the longitudinal momentum fractions of the quark and antiquark given by a. That 
reflects the usual factorisation of transverse and longitudinal degrees of freedom in high 
energy reactions. But also for this assumption a general proof is not known. 

Note that we have assumed here a dependence of the reduced matrix element on 
the squared energy s rather than on Bjorken-x. We will explain our choice and discuss 
it in more detail in section 6.4 below. 

With the assumptions (iii)-(v) we find for the spin averaged T-matrix element at 
high energy 



|^5,.,(D®(q,a,RT,A',A),p(p,5')|r|D('?)(q,a,RT,A',A),p(p,5)) 



2 



where T^^l is the reduced T-matrix element and involves only contributions coming 

from T'"-*. The fact that 7^^^ is a function of R^ and s only is an immediate conse- 
quence of assumption (v). The more conventional definition of the reduced T-matrix 
element would include all factors depending on the continuous parameters of the T- 
matrix element. Here we have explicitly taken out the two singular factors given by the 
two delta-functions involving a and R^^ in order to have a reduced T-matrix element 
that can as usual be assumed to be a smooth function of its arguments. Alterna- 
tively, we could have chosen to approximate the delta-functions involving a and Rj^ in 
(120) by strongly peaked but smooth functions and to absorb them into the definition 
of the reduced matrix element. That would correspond to requiring an approximate 
diagonality in assumption (iii). 

We can now obtain the reduced cross section for the scattering of a dipole state on 
an unpolarised proton from the reduced T-matrix element via the optical theorem, 

ai^) (R|, s) = \ Imr,(^i(R2., .) > . (121) 

Recall that we have assumed that quarks have a mass-shell in assumption (i) . Therefore 
we can now conclude that this reduced cross section is non-negative since it can due to 
that assumption be related to a physical three-to-three scattering process qqp -^ qqp. 
In this sense we can further below interpret the reduced cross section 0"^^^ as the cross 
section for a scattering of a dipole on an unpolarised proton. Note that the reduced 
cross section can naturally be expected to depend on the quark flavour q, for instance 
via the quark mass niq. 

Instead of considering the original dipole states D we could also consider the cor- 
responding T-matrix element of smeared dipole states D as introduced in (103). Then 
we obtain from (120) for large s 

s',s 

= ^mh'yhx^^MT,^) (122) 



with 



f}^{K'T, s) = J SRt |/t(Rt - Rt)|' r}^{^l, s) . (123) 
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We can interpret 7^^^ as the genuine proton spin averaged T-matrix element for smeared 
dipole states and can assume that it is a smooth function of its arguments. Note that 
in the matrix element (122) the incoming and outgoing dipole states need to have the 
same a and Ry in order to arrive at (123). We can then again invoke the optical 
theorem to obtain the total cross section for a smeared dipole state, a hadron analogue, 
scattering on an unpolarised proton for large s as 



--^m^s) 



1 



-(<?), 



2 M , 



Imr^^^,^(R^, .) = / d'RT |/r(RT - Rt)|^ a'^^^^ ^) • (124) 



For a narrow smearing function /j^ we get o'j:^^(Ry, s) « ^rcd(-'^T' •*)' '^^^ is (121) gives 
a correctly normalised hadron-analogue cross section. 

Returning now to the unsmeared dipole states we obtain with (120) from (118) 
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This gives us for the hadronic tensor of (110) at high energy 
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Accordingly, we find with (111) and (112) for o"t and a^ at high energy 
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where we have defined 
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(129) 
(130) 



We recall that in (130) we have to use the polarisation vector e^ (56) since only then 
do we get the correct asymptotic expression, as was explained in section 4. Here and 
in the rest of the paper we only deal with the asymptotic expressions and no longer 
indicate this by a subscript. Note that wY^ and w)^ depend only on Rt but not on 
the orientation of Ry due to (54) and (61). 
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With (127)- (130) we have finahy arrived at the standard formulae for the dipole 
picture used extensively in the literature, see for example [17]. We have tried to spell out 
all the assumptions that enter and that should be carefully examined when one wants to 
draw conclusions, for instance concerning saturation, from a supposed 'unitarity limit' 
for a^cdi^T^ ^)- ^^ consider it an important task for the future to test the assumptions 
and to quantify their accuracy. 

Note that in the dipole formulae (127) and (128) Rt appears formally as an inte- 
gration variable. But the interpretation of these formulae as a factorisation into the 
square of a photon wave function and a cross section for the scattering of a dipole of 
transverse size Rt on a proton is very natural in view of the relation of the latter to the 
corresponding T-matrix element, see (120) and (121). We will therefore simply call a^'^j^ 
a dipole cross section in the following. Due to the same relation to a T-matrix element 
involving a particular quark flavour it is also natural that this dipole cross section still 
depends on that quark flavour q. 

Next we discuss two issues that we have postponed earlier in this section. 

6.3 Proton spin dependence of the dipole cross section 

The first point is the dependence on the proton spin. So far we have considered only 
the scattering of transversely or longitudinally polarised photons on an unpolarised 
proton, and we have obtained the usual formulae of the dipole picture for the total 
cross sections of these processes. As a consequence, we needed to make our assumption 
(iii) only for the proton spin averaged T-matrix element (119) rather than for the 
actual T-matrix element which occurs in (118). Nevertheless we have formulated our 
assumption (iii) for the latter, which amounts to making a stronger assumption. It 
now seems very interesting to think about possible tests of the dependence of that 
matrix element on the proton spin. To our knowledge, all models proposed so far for 
the dipole cross section a^^d have been used only for calculating scattering processes 
on unpolarised protons. Some of these models are based on an additional assumption 
about the proton, namely that it can be considered at high energies as a superposition 
of colour dipoles described by a suitable weight function for their transverse size and 
orientation. The scattering process can then be expressed in terms of the cross section 
of two colour dipoles. In phenomenological applications one hence only needs to make 
a model for that simpler cross section. For a typical approach along these lines see for 
example [18]. We emphasise that this assumption goes beyond the dipole picture as we 
have derived it here, and it is in our opinion far from obvious that the proton can be 
approximated by a collection of dipoles. It would of course be very interesting to see 
if that assumption can be justified using the techniques developed in the present work. 
Let us now for a moment adopt that additional assumption. One expects the dipole- 
dipole cross section and the corresponding T-matrix element to be symmetric under 
the exchange of the two dipoles at least in the large- A^c limit in which a colour dipole 
can be viewed as a quark-antiquark pair. In assumption (iii) we have assumed that 
the T-matrix element in (118) is proportional to the unit matrix in the space of spin 
orientations of the quark and antiquark. Accordingly, one would now have to assume 
that the dipole-dipole T-matrix element is also proportional to the unit matrix in the 
space of spin orientations corresponding to the dipole in the proton. As an immediate 
consequence one finds that in such a picture the polarised proton structure function 
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gi vanishes. Since this is in contradiction with the experimental findings (for a review 
see [19]) one might conclude that at least the additional assumption does not hold, but 
that would be premature. On the contrary, that result would be in agreement with the 
general finding that the polarised structure function gi is suppressed at high energy with 
respect to the unpolarised structure function F2. It has been established for example in 
Regge theory that gi is suppressed with respect to F2 by a power of the energy [20], and 
a similar picture emerges also in perturbation theory both in extrapolations of Altarelli- 
Parisi evolution [21] to small Bjorken-x [22, 23] and in calculations taking into account 
double-logarithmic contributions beyond Altarelli-Parisi evolution [24, 25]. In many 
steps of our arguments leading to the dipole picture subleading terms at high energy 
have been neglected. Thus it is consistent that the leading level result for an observable 
that starts at a subleading level is zero. Turning this argument around the polarised 
structure function gi might offer an opportunity to study the structure and the size of 
subleading terms at high energies. The most important question is of course whether it 
is possible to find a consistent description of gi (or of any other observable that starts 
at a subleading level in the expansion in inverse powers of the energy) in the dipole 
picture. This problem could be studied independently of the additional assumption 
that the proton can be considered as a superposition of colour dipoles. (Note, however, 
that that additional assumption is at least consistent with the known results about ^fi.) 
Further studies along these lines might either show that also subleading terms can be 
factorised according to the dipole picture, or tell us about a potential breakdown of the 
dipole picture for subleading terms. In the latter case one could even hope to quantify 
the size of potential corrections to the dipole picture in general. 

6.4 The energy variable of the dipole cross section 

The other important issue that we have not yet discussed in detail is the question 
whether the dipole cross section cji-cd should depend on the squared energy s or on 
Bjorken-x in addition to the dependence on the dipole size Rt- Since at high energies 
XBj = Q'^/s the latter possibility would imply a dependence of the dipole cross section 
on Q^ while such a dependence would be absent if the former possibility were chosen. 
Phenomenological models for the dipole cross section have been proposed based on both 
possibilities, and both can lead to a satisfactory description of the available data. Over- 
all, models based on an XBj-dependent dipole cross section appear to be more popular 
at present. A model with an s-dependent dipole cross section has been constructed for 
instance in [26] while a prominent example for an XBj-dependent dipole cross section 
is the model of [27]. In our approach we have obtained the dipole cross section from a 
T-matrix element for the scattering of a dipole state on a proton, see (120) and (121). 
But as we have discussed in section 5.2 the dipole states entering this matrix element 
are independent of the photon virtuality Q^ at high energy as can be seen from the 
expression for their invariant mass (99), see also the corresponding formulae (105) and 
(106) for smeared dipole states. The relevant T-matrix element thus corresponds to a 
scattering process in which the initial and final state is independent of Q^. According 
to this observation we have chosen our dipole cross section a^^d to depend only on s 
rather than on xbj in our assumption (v). The main motivation for choosing an xbj- 
dependence in many phenomenological models comes from studying the simultaneous 
limit of large s and large Q^. In this limit the cross section for photon-proton scattering 
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should agree with the perturbative result obtained in the double-leading logarithmic 
approximation (DLLA) [28]. As was discussed already in [29] the dipole cross section 
can in that limit be identified up to factors with the gluon density of the proton. The 
latter naturally depends on xbj and Q^ and therefore suggests to choose the dependence 
on xbj also at lower Q^ . Can we decide which variable is the correct one for the dipole 
cross section? 

Let us first recall that Ry in (127) and (128) is an integration variable. We can 
therefore perform a change of variables 

Rt = ^R^, (131) 

where Qo is some fixed momentum scale. Expressing the integrands of (127) and (128) 
in terms of the new integration variable R^ we obtain formulae similar to the dipole 
model that contain a factor 

-S(R^^Q,-BJ)=ai1(f R?>-) (132) 

which we could define as a reduced cross section that depends on R^ , Q and xbj instead 
of Ry and s. However, in this case R^ is no longer the physical size parameter of the 
colour dipole. Moreover, the substitution (131) would affect also the other factor in 
the integrands of (127) and (128) and hence spoil the interpretation of the integrand as 
a product of an actual dipole-proton cross section with the square of the perturbative 
photon wave function. We therefore discard the mathematical possibility of changing 
the integration variable from Rt to R^. 

A discussion of s versus xbj dependence of the dipole cross section was given in chap- 
ter 9 of [17] and our remarks below follow the same lines. For large photon virtualities 
Q^ one finds that the integrals in (127) and (128) receive the dominant contribution 
from small R^ due to the shape of the weight factors w^ i^- For small Rt the dipole- 
proton cross section should become small due to colour transparency, and the typical 
behaviour to be expected at small Rt is (Tred(R'T) s) oc i?|,. Hence a typical ansatz for 
small Rt is 

^S(^?^,^) = ^?^/(^), (133) 

where / is a dimensionless scaling function. Naive scaling arguments suggest this and 
in QCD we expect (133) to hold up to logarithmic corrections in Rt- Inserting (133) in 
(127) and (128) we find that the explicit factor J^ together with the weight functions 

w^ j^{Rt,Q'^) and an additional factor Rt from the integral measure cPRt leads to a 
pronounced maximum in the integrands situated at 

^Tmax=Q2+4^2 (13^) 

with C^ ~ 4, see p. 273 of [17] and figure 9 in section 6.5 below. That is, for large Q^ 
mostly dipoles of size Rt max contribute. We may therefore make the replacement 



^f(^)^-illiR'T,XB,) (135) 
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for Q^ ^ Ami. Hence it is only due to the behaviour of the other factors in the 
integrands and the exphcit factor i?|, in (133) that an s-dependent a^^A and an xbj- 
dependent iTred lead to the same result at large Q^ . It is therefore not possible to 
establish an XBj-dependence of the dipole cross section in general based on this argu- 
ment. As we have described above, in our approach an s-dependence occurs naturally. 
We cannot exclude that a different derivation of the dipole picture based on a different 
set of assumptions is possible in which an XBj-dependence occurs. If one wants to main- 
tain the interpretation of iTred as an actual cross section involving a quark-antiquark 
colour dipole state, however, one would need to define the dipole states in such a way 
that they depend on Q^ in leading order in the expansion in inverse powers of s. This 
would in particular exclude the natural definition (95) which we have used. 

In this context we should add a further remark. It would of course be interesting 
to establish the relation of our derivation of the dipole picture with the perturbative 
picture of the DLLA step by step. But there is a potential problem. In the present 
paper we study the limit 

s 

whereas in the DLLA one takes 



oo , Q fixed , (136) 



s ^ cxD , Q^ ^ cxD , — fixed. (137) 

s 

It is not at all clear and not necessary that first taking the limit (136) and then letting 
Q^ become large will lead to the same result as taking the limit (137). It remains to 
be seen whether this makes the comparison of our derivation with the perturbative 
approach intrinsically difficult. 

6.5 Density for the photon wave function in leading order 

In the dipole picture the photon-proton cross section is obtained as the convolution of 
the square of the photon wave function with the dipole cross section. The latter cannot 
be calculated from first principles at present and needs to be described by models in 
phenomenological applications of the dipole picture. The photon wave function, on 
the other hand, is known explicitly in perturbation theory. We therefore find it use- 
ful to discuss here some known aspects of the photon wave function and to illustrate 
some of its properties. In the next section we will then derive some phenomenologi- 
cal consequences of the dipole picture which will be based only on the properties of 
the perturbative photon wave function but will not depend on any particular model 
assumptions about the dipole cross section. 
In the following we consider the densities 

4^)(a,i?T,g') = 5^|vi';lv(a,RT,Q)f (138) 

A,A' 
= ^ «cmQ^ { {o? + (1 - af\ el{K^(e^RT)f + m\{K^{e^RT)f] 
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and 



v'^l\a,RT,Q^) 



^ v;Xa'(«>Rt,q) 



A,A' 



TT^ 



«emQ'Qla(l - a)Y[Ko{eqRT)y 



(139) 



for transversely and longitudinally polarised photons as obtained from the respective 
photon wave functions in leading order, see (54), (61) and (50). Upon integration over 
the longitudinal momentum fraction a these densities give the functions wt and wl in 
(129) and (130) that occur in the dipole formulae (127) and (128), 



w!j!J^{Rt, Q^) = da 4% (a, Rt, Q^) ■ 



(140) 



For our numerical results we assume the light quarks (q = u, d, s) to be massless while 
for the heavy quark masses we use rric = 1-3 GeV and nib = 4.6 GeV. 

We start with transversely polarised photons. Figure 5 shows the dependence of the 
density w^ /(acmQg) on the longitudinal momentum fraction a for fixed Rt = 1 GeV~ 
and for three different values of the photon virtuality Q^ = 1, 10, 100 GeV . The left 
panel is for massless {u, d, s) quarks while the right panel is for c- and fe-quarks. Note 
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Figure 5: Dependence of the density 



Tm^T 



vi^Ua,RT,Q^) 



of transversely polarised 



photons (see (138)) on a for Rt = 1 GeV~^ and for three different values Q^ = 
1, 10, 100 GeV^, plotted in units of GeV^; on the left for massless {q = u, d, s) quarks, 
on the right for c-quarks (upper lines) and 6-quarks (lower lines). 

that for transversely polarised photons there is a nonvanishing contribution from the 
end points a = and a = 1 even for very large photon virtualities Q^. As can be 
seen also from (138) the value of Vj! at these end points depends on the mass of the 
quark flavour under consideration. The density has a minimum at a = 1/2 for all 
values of Q^ and the relative size of the density at this minimum compared to the end 
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points shrinks with increasing Q^ . That means that a symmetric distribution of the 
longitudinal momentum of the photon between the quark and antiquark is the least 
likely configuration, and this effect becomes more pronounced for larger Q^ . The size of 
the density decreases rapidly with increasing Q^ for all a with the exception of the end 



points. Figure 6 shows the dependence of v)^ /{a, 

Q^ = 10 GeV^ and for three different values of the dipole size Rt 



Q'l) on Q for fixed photon virtuality 



0.1,l,5GeV-\ 
again for massless quarks on the left and for massive quarks on the right. While for 
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Figure 6: Dependence of the density 



:m^T 



v^^ha,RT,Q^) 



of transversely polarised 



photons (see (138)) on a for Q^ = lOGeV and for three different values Rt = 
0.1, 1, 5GeV~^, plotted in units of GeV^; on the left for massless {q = u,d,s) quarks, 
on the right for c-quarks (upper lines for each line type) and 6-quarks (lower lines for 
each line type). 

Rt = 0.1 GeV^ the densities for c- and 6-quarks cannot be distinguished when plotted 
on a logarithmic scale like in figure 6 their sizes become very different for larger Rt- 
Already for Rt = 5 GeV^ the density for the 6-quark is so small that it is far below 
the range shown in the figure. For all quark flavours the dependence of the density on 
a changes with Rt in a similar way as it changes with Q, in particular it decreases 
rapidly with increasing Rt for all momentum fractions a. For massless quarks this 
similarity can be understood immediately based on dimensional arguments. Namely, 
for massless quarks the dependence of the density on Q^ is related to the dependence 
on Rt because the dimensionless quantity R^v^ can depend only on the product QRt 
by dimensions. However, we have plotted the dimensionful density vt in our figures 
and hence the values of the density at the end points a = and a = 1 are different for 
different values of Rt- 

Next we turn to longitudinally polarised photons. In figure 7 we show the de- 
pendence of the density f]J' /(oernQo) o^ ^^^ longitudinal momentum fraction a for 
fixed Rt = 1 GeV~ and for three different values of the photon virtuality Q^ = 
1, 10, lOOGeV^. The left panel is again for massless {u,d,s) quarks while the right 



36 



panel is for c- and 6-quarks. In the case of longitudinally polarised photons the density 
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Figure 7: Dependence of the density 



aQq 



y^ {a, Rt, Q^) of longitudinally polarised 

photons (see (139)) on a for Rt = 1 GeV~ and for three different values Q^ = 
1, 10, 100 GeV^, plotted in units of GeV^; on the left for massless {q = u, d, s) quarks, 
on the right for c-quarks (upper lines) and 5-quarks (lower lines). 

vanishes at the end points a = and a = 1 for all photon virtualities Q^ and for all 
quark flavours. For small photon virtualities the density has a maximum at a = 1/2. 
For larger Q^ the symmetric distribution of the longitudinal momentum of the photon 
between the quark and the antiquark becomes less likely and there appears a local 
minimum in the density at a = 1/2. At the same time two maxima of the density 
develop close to the end points a = and a = 1. For massless quarks the overall size 
of the density decreases rapidly with increasing Q'^, except for small regions close to 
the end points where the maxima appear at larger Q^. For massive quarks that overall 
decrease with Q^ sets in only for Q^ > rni. Figure 8 shows the dependence of the 
density on a for fixed photon virtuality Q^ = lOGeV and for three different values 
Rt = 0.1, 1, 5 GeV~ , again for massless quarks on the left and for massive quarks on 
the right. The size of the density decreases with increasing Rt, and for Rt = 5 GeV~ 
the density for 6-quarks is already so small that it is far below the range shown in the 
figure. For massless quarks the dependence of the density on Rt can again be deduced 
from the dependence on Q^ via the dimensionless quantity Rt^l which depends only 
on QRt- For large Rt the dependence of the density for massive quarks on Rt follows 
a pattern similar to that of massless quarks. 

Finally, we consider the dimensionless quantity QR^w^\^ for transversely and lon- 
gitudinally polarised photons which occurred in our discussion of the energy variable 
of the dipole cross section in section 6.4. There we pointed out that the densities w)^^ 
occur in the dipole formulae (127) and (128) together with a factor Rt from the integral 
measure and another factor of R^ coming from the dipole cross section at small Rt, see 

(133). Multiplying by Q leads us to the dimensionless quantity QRj^Wr^]^. We recall 
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Figure 8: Dependence of the density 



,(9) 



aQq 



a, Rt, Q"^) of longitudinally polarised 



photons (see (139)) on a for Q^ = lOGeV and for three different values Rt = 
0.1, 1, 5GeV~^, plotted in units of GeV^; on the left for massless [q = u,d,s) quarks, 
on the right for c-quarks (upper lines for each line type) and 6-quarks (lower lines for 
each line type). 



in (138) and (139) by integrating over a, see (140). 

(9) 



that ti'j?^ is obtained from vj!^ 

Figure 9 shows the dependence of QRj^w^'j^ on Rt for massless quarks {q = u,d,s 
and for three different photon virtualities Q^ = 1, 10, 100 GeV . The left panel is for 
transversely polarised photons and the right panel for longitudinally polarised photons. 
The figure shows that there are indeed pronounced maxima, and their position i?Tmax 
scales with l/Q, see also (134). For transversely polarised photons the maxima turn 
out to be wider than for longitudinally polarised photons. Since for massless quarks 
the dimensionless QRj^w^ j^ can depend only on QRt the value of this quantity at the 
maximum is the same for all Q^ here. 

7 A bound on R = (Tl/ctt from the dipole picture ^ 

In the previous section we have obtained the dipole picture of high energy scattering 
based on a number of assumptions. It is obviously important to test these assumptions 
and to find the range of kinematical parameters in which they hold, or in other words, 
the range in which the dipole picture can be applied. 

In the dipole picture the cross section for photon-proton scattering is written as a 
convolution of the square of the perturbative photon wave function with the dipole- 
proton cross section, see (127) and (128). The former can be determined explicitly 
in perturbation theory (see section 4) while the latter cannot be calculated from first 
principles at present. Especially in the interesting region of intermediate and low photon 
virtualities one therefore has to retreat to models for the dipole cross section which can 



iSome of the results in this section have been presented in abridged form also in the Letter [30]. 
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Figure 9: Dependence of the diniensionless quantity QRj^wl^^ on Rt for niassless 
{q = u,d,s) quarks for three different values Q^ = 1,10, 100 GeV ; on the left for 
transversely polarised photons, on the right for longitudinally polarised photons. 



then be tested against the available data. Although the required agreement with the 
data places constraints on the possible models for the dipole cross section there is still 
considerable freedom in these models. Based on such models it is therefore difficult to 
test the assumptions on which the dipole picture is based and to explore their potential 
limits. It would hence appear favourable to find observables which can be used to test 
the assumptions of the dipole picture in a way which does not depend on any particular 
model for the dipole cross section. 

In the present section we want to study such an observable, namely the ratio R of 
the cross sections for longitudinally and transversely polarised photons, 



Ris,Q^ 






(141) 



In the following we will derive a bound on R from the dipole picture which will involve 
only the properties of the perturbative wave functions of longitudinally and transversely 
polarised photons. 

First we recall that the quantities Wrp and wj^ in (129) and (130) are positive since 
they are obtained as integrals over the squared modulus of the perturbative photon wave 
function. Further we recall that the dipole cross section is non-negative according to 
our assumptions made in section 6.2, see (121). We then start from the obvious relation 



w 



(<?) 



iRT,Q^ 



W^ 



(<?) 



< max 



{Rt,Q') «'^^ 



w^^\Rt,Q') 



(142) 



w. 



where the maximum is taken over all dipole sizes Rt and over all quark flavours q. We 
then multiply both sides of this relation by the positive quantity w!^ {Rt, Q"^) and by 



the dipole cross section (J^J^{R' 



T' ' 



which is non-negative. Next we integrate over Ry 
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and sum over all quark flavours which gives with the dipole formulae (127) and (128) 



o'l{s,Q ) < crT{s,Q ) niax 



,(<?) 



9.«T w)§\Rt,Q'^) 



(143) 



An analogous relation is readily obtained for the minimum of wy^ /wj! , such that we 
have 



w 



(?) 



mm 

q,RT 



{Rt,Q^ 



w, 



(<?), 



< R{s,Q^) < max 



w 



(?) 



iRT,Q^ 



(<?), 



(144) 
.^\Rt,Q^) ' ' ■ «'^^ w}^\Rt,Q^) 

With this relation we have obtained an upper and a lower bound on the ratio R{s, Q'^) 
which should be respected for all s and Q^ for which the dipole picture is valid or ~ 
phrased differently - for which the assumptions and approximations (i)-(v) hold which 
we have discussed in section 6.2. Note that the bounds (144) depend only on the photon 
virtuality Q^ but do not depend on s. This is of course expected since the bounds 
involve only the photon wave function but not the dipole cross section cJi-cd(-Rr' '^)' ^^^ 
only the latter contains the energy dependence of cjt and ai according to the dipole 
formulae (127) and (128). 

Similarly, we can obtain bounds for the ratio of the cross sections for the production 
of a heavy quark flavour, q = c or q = b, in deep inelastic scattering. 



mm —if- < Rq{s,Q ) 



M)i 



,(i)i 



— < max 



4'^(s,Q2 



Si) 



Rt w}^>{Rt,Q^) 



(145) 



Let us now determine the numerical values of the bounds derived here. In figure 
10 we show the ratio w^ {Rt,Q^)/wj^ [RtjQ"^) for massless {q = u,d,s) quarks as 

= 10 GeV . This ratio has a maximum at 




,(<?)( 



,(«) 



Figure 10: The ratio w"^' {Rt,Q^)/w^j!' {Rt,Q^) as a function of Rt for Q^ = 10GeV^ 
for massless {q = u, d, s) quarks. 
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which its value is 0.37248, and asymptotically it approaches zero. Since w^ /wj! is a 
dimensionless quantity it can for massless quarks depend only on QRt such that also for 
different photon virtualities Q^ we find a maximum with the same value. For massive 
quarks {q = c, b) the ratio w^ /w^ as a function of Rt also exhibits a maximum with 
a value which now depends on Q^, but for all Q^ we find that this value is smaller than 
for massless quarks. 

We can therefore conclude that the upper bound on R in (144) gives 



R{s,Q'^) < 0.37248, 



(146) 



while the lower bound on R in (144) reduces to the trivial statement that R is non- 
negative. Also the lower bound on the ratios Re and Rb in (145) is trivial. The upper 
bounds on Re and Rb resulting from (145), on the other hand, are nontrivial and depend 
on Q^. This dependence is shown in figure 11 together with the upper bound (146) on 
R. 
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Figure 11: The upper bound on R = ai^jOT resulting from the dipole model as a 
function of Q^ . The solid line is for massless (u, d, s) quarks and constitutes the 



bound (144). The lower lines are upper bounds on R^ 
respectively, see (145). 
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^L I^T ^^^ Rb = <^L /c. 
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In many practical applications of the dipole picture one chooses to introduce phe- 
nomenological masses for the light [q = u, d, s) quarks. In [27] for example a light quark 
mass of mq = 140 MeV is chosen, and this mass is implemented in the perturbative 
photon wave function. Such a modified photon wave function results in a corresponding 
modification of the bound on R which then becomes Q^-dependent. For rUg = 140 MeV 
this bound is shown on a linear scale in figure 14 below. With the logarithmic Q^-scale 
in figure 11 it would correspond to a curve of the same shape as the bound on Re or 
Rb but shifted to the left such that its value at Q^ = 1 GeV would be 0.318. 

A popular extension of the dipole picture as we have discussed it so far is to diffrac- 
tive deep inelastic scattering. In this extension one usually writes for the diffractive 
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photon-proton cross section [31] 

{s,Q') = -^Y.Jd'RTW^^j^{RT,Q') [a'^l\{Kl,s)y , (147) 



daff. 



dt 



t=o q 



where the dipole cross section cired is the same as in the usual dipole picture. The dipole 
formula for diffractive scattering (147) therefore relates diffractive to inclusive photon- 
proton scattering. Accepting (147) as a valid description of the diffractive cross section 
one finds that exactly the same bounds as obtained above apply also to diffractive 
photon-proton scattering, both for the total diffractive cross section and for diffractive 
production of heavy quarks. In particular, we have 



Zl 

dt 



t=o 



-rdiff 
dt 



t=o 



< 0.37248 . (148) 



We emphasise, however, that the extension of the dipole picture to diffractive scat- 
tering (147) goes beyond what we have discussed so far in the present work. It is in 
particular not clear and in fact not necessary that the dipole formula (147) for diffrac- 
tive scattering follows from the same assumptions (i)-(v) which we have used to obtain 
the dipole picture for inclusive deep inelastic scattering. It is an open question whether 
it is possible to establish the dipole picture of diffractive photon-proton scattering based 
on the techniques presented here. 

The dipole picture has been applied also to a number of exclusive reactions, a typical 
example is diffractive vector meson production. Here the wave function of the outgoing 
photon is replaced by the wave function of the produced vector meson, and the latter 
usually involves some model assumptions. We should point out that our bounds on 
R cannot be applied directly to such processes. This is because our derivation of the 
bounds relies on the positivity of both factors under the integrals in the dipole formulae 
(127) and (128). In particular, it relies on the positivity of the square of the photon 
wave function which occurs due to the incoming and the outgoing photon. In other 
reactions that factor is changed and positivity is in general not guaranteed. 

We now return to inclusive deep inelastic photon-proton scattering and discuss the 
experimental data on R. The presently available data on R at high energies have been 
obtained by the NMC [32], CCFR [33], E143 [34], EMC [35] and CDHSW [36] collab- 
orations. The scattering processes used to extract R include not only e p scattering 
but also a variety of other scattering processes including muon and neutrino scattering 
on nuclear targets. In some of these processes one might in general expect additional 
caveats concerning the applicability of the dipole picture. For lack of further experimen- 
tal data we nevertheless include the corresponding data points here. At high energies 
R is expected to be independent of the target, as is also confirmed by the data within 
the experimental errors. For further details we refer the reader to the experimental 
publications. 

We want to consider only data points at sufficiently high energy which we choose 
to be data with xbj < 0.05. Figure 12 shows all available data in this energy region 
together with the bound (146) resulting from the dipole picture. The few data points 
available at xbj < 0.01 are shown as full points while those with 0.01 < xbj < 0.05 are 
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Figure 12: Comparison of experimental data for R = ui^jaT in the region xbj < 0.05 
with the bound (146) resulting from the dipole picture. Full points correspond to data 
with a;Bj < 0.01, open points are data with 0.01 < xbj < 0.05. 

shown as open points. The data have rather large errors, but by and large we can say 
that they respect the bound resulting from the dipole picture. For Q^ below 2 GeV , 
however, there appears to be the tendency that the data come close to the bound. Note 
the interesting fact that data very close to the bound (146) could be accommodated 
in the dipole picture only if the dipole cross section a^^A were strongly peaked in Rt 
around the maximum of w]^ /w)^ - which would appear to be very unlikely in reality. 
Hence already data close to the bound could be interpreted as an indication of the 
breakdown of the dipole picture. The errors of presently available data are far too 
large to draw any firm conclusion about the assumptions underlying the dipole picture 
here, but we might take them as an indication for the values of Q^ below which we 
should be careful in interpreting scattering processes in terms of the dipole picture only. 

Unfortunately, no direct measurements of F^ and R have been performed at HERA 
so far. A discussion of the available indirect determinations of Fl in view of our 
bound will be given elsewhere. Direct measurements of Fl at HERA are planned and 
will hopefully lead to a better determination of i?, allowing for a stringent test of the 
validity of the dipole picture at low Q^ . Also a potential future upgrade of RHIC to an 
electron-ion collider eRHIC would offer the possibility to obtain a clear picture of R at 
high energies and low Q^, and thus to establish a conclusive test of the dipole picture 
in this region. 

Let us finally see which behaviour of R{s,Q'^) is predicted in the framework of a 
typical model for the dipole cross section a^cd- As a prominent example we choose the 
Golec-Biernat-Wiisthoff model [27]. In that model the dipole cross section is assumed 
to be a function of Rj<, xbj and Q^ and not only of i?j. and s. This might pose a problem 
in the context of our discussion in section 6.4 but for the purpose of illustrating the 
behaviour of i? in a typical dipole model we can disregard this issue. We should point 
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out in this context that our bound (146) is completely independent of the variables 
on which the dipole cross section depends. In the Golec-Biernat-Wiisthoff model the 
dipole cross section Ured is given by 






with 



'-^"p(-5| 



A 

-2 



(149) 



i?^(x) = I — I GeV"^ , (150) 

where the three parameters of the model are chosen to be ctq = 23 mb, A = 0.29, 
xo = 3 • 10"'^. The variable x in (150) is obtained from xbj via 

x = x^Al + -^\ , (151) 

due to which the dependence of o"gbw on the quark flavour and on Q^ enters. We use 
again rric = 1.3 GeV for the charm quark mass. In order to study the photoproduc- 
tion limit Q^ ^ a phenomenological mass of m^ = 140 MeV is chosen for the 
light quarks {q = u, d, s) which serves as a regulator for logarithmic divergences that 
would otherwise occur in this limit. In the Golec-Biernat-Wiisthoff model [27] that 
phenomenological light quark mass is implemented not only in the dipole cross section 
via (151) but also in the photon wave function. Note, however, that after inclusion of 
the effects of DGLAP evolution a better fit to the data is obtained with a vanishing 
light quark mass [37]. 

In order to compare with our bound (146) we first insert the dipole cross section 
^GBW (149)-(151) with the light quark mass m'^^^ = 140 MeV into the dipole formulae 
(127) and (128) but keep the light quarks massless in the photon wave functions. The 
resulting behaviour of i? as a function of Q^ is shown in figure 13 for four different values 
of XBj, xb] = 0.05,0.01, 10^'^, 10~^, together with the bound (146). The prediction for 
R depends only weakly on xbj and stays significantly below the bound. 

Next we consider the original Golec-Biernat-Wiisthoff model, that is we implement 
the phenomenological light quark mass m'^^^ = 140 MeV not only in the dipole cross 
section but also in the photon wave function. The resulting Q^-dependence of R for 
XBj = 0.01 is shown as the solid line in figure 14. The dashed line in that figure is the 
one which we obtained before in figure 13 by using the light quark mass m^^^ only 
in the dipole cross section. At low Q^ the phenomenological quark mass leads to a 
significant suppression of R while for Q^ > 3 GeV the light quark mass in the photon 
wave function does not have a significant effect. A comparison with figure 12 shows that, 
given the large errors of the data, both curves are in agreement with the measurements 
of R. The dot-dashed line in figure 14 is the bound obtained from a photon wave 
function with a light quark mass of 140 MeV as we discussed it already earlier in this 
section. For photon virtualities Q^ below about 2 GeV the bound is significantly 
below the original bound (146) which was obtained for massless light quarks and which 
is shown in the figure as the dotted line. 

More precise data on R might in the future serve as an additional test of the 
quality of various models for the dipole cross section. It would therefore be interesting 
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Figure 13: The ratio R = cfl/'^t as obtained from the Golec-Biernat-Wiisthoff model as 
given by (149)-(151) but with vanishing hght quark mass in the photon wave function. 
The sohd hues are for four different values xbj = 0.05, 0.01, 10~^, 10~^ from bottom to 
top at large Q^. The dashed line is the bound (146) derived from the dipole picture. 
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Figure 14: The ratio R = cfl/'^t as obtained from the Golec-Biernat-Wiisthoff model 
in comparison with the bound on R. The solid line results from the model (149)-(151) 
for XBj = 0.01 with the light quark mass 140 MeV in the dipole cross section and in the 
photon wave function, the dashed line is obtained if that light quark mass is used only 
in the dipole cross section. The dotted line is the bound (146) and the dot-dashed line 
is the modified bound as obtained for a light quark mass of 140 MeV. 
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to study in more detail the behaviour of R also in other models for the dipole cross 
section. An interesting aspect of such studies could be in particular the choice of a 
suitable phenomenological light quark mass which is required for the photoproduction 
limit or is motivated by the expected breakdown of the perturbative description of the 
photon wave function at large dipole sizes. For a consistent description one should 
clearly use the same light quark mass in the photon wave function and in the model 
for the dipole cross section. Since the bound on R depends on the mass of the light 
quarks in the photon wave function especially at low Q^ it would even be possible to 
use more precise data on R in that region to obtain bounds on possible choices of a 
phenomenological mass for light quarks in the photon wave function and hence also in 
models for the dipole cross section. 

8 Conclusions and outlook 

In this paper and in the companion paper I we have studied how the dipole picture of 
high energy scattering can be derived in the framework of a genuinely nonperturbative 
formulation of photon-nucleon scattering. We have analysed in detail the Compton 
amplitude for real and virtual photons and have isolated the contributions to it which 
are leading at high energies. We have taken into account properly the renormalisation 
of these contributions and have studied their gauge invariance emphasising in particular 
its relevance for the definition of the perturbative photon wave function. We have then 
identified the approximations and assumptions which are necessary to arrive at the 
usual dipole picture. 

At high energies we obtain indeed the dipole picture, but also find two important 
additional contributions. One of them corresponds to rescattering effects of the quark 
and antiquark forming the colour dipole. The other additional contribution corresponds 
to diagrams in which the incoming photon and the outgoing photon do not couple to the 
same quark loop. Both of these terms are of higher order in the coupling constant a^ 
but are not suppressed by powers of the energy. In lowest order in perturbation theory 
these two additional terms do not contribute, and one is left with only the usual dipole 
picture. At low photon virtualities on the other hand we expect that both additional 
contributions can have sizable effects. We expect for example that the rescattering term 
will play an important role in the transition to low photon virtualities for which vector 
meson dominance applies. The rescattering term should in fact generate the relevant 
diagrams for the formation of quark-antiquark bound states. It would therefore be 
very interesting to see whether our results can be used to obtain a better picture of 
the transition from a perturbative photon wave function at high virtualities to vector 
meson dominance at low virtualities. 

A crucial question in all phenomenological applications of the dipole picture con- 
cerns the range of its applicability, or in other words the kinematical region in which 
the approximations and assumptions leading to the dipole model hold. We have derived 
a bound on the ratio R of the cross sections for the scattering of longitudinally and 
transversely polarised photons off hadrons which is relevant in this context. The bound 
is obtained only from the perturbative photon wave function and is independent of any 
particular model for the dipole-proton cross section. Therefore it allows for a model 
independent test of the dipole picture. Unfortunately the presently available data are 
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not precise enough to restrict the range of apphcabihty of the dipole picture. However, 
the data appear to come close to the bound below photon virtualities of about 2 GeV , 
suggesting that at least in that region one should be careful in interpreting high energy 
scattering data in terms of the dipole picture only. In [30] we have presented further 
bounds on ratios of nucleon structure functions F2 taken at the same energy y^ but 
at different photon virtualities Q^ resulting from the dipole picture. They are suitable 
for further constraining the range of applicability of the dipole picture. 

The dipole picture is often discussed in the context of the space-time picture of 
high energy photon-hadron scattering in the target rest frame. Here the photon splits 
into a long-lived quark-antiquark pair which then interacts with the hadron. Our 
study confirms this space-time picture and supplements it with the renormalisation 
of the photon-quark-antiquark vertex. The space-time picture is closely related to 
the notion of the light-cone wave function of the photon, suggesting a Hamiltonian 
picture of the evolution of the photon wave function and of the scattering process. We 
emphasise that we have obtained the dipole picture in a purely Lagrangian approach 
here. This allowed us in particular to implement the renormalisation procedure in a 
simple and transparent way. In a Hamiltonian approach a natural extension of the 
dipole picture is to include higher Fock states in the photon wave function. In a first 
step one would include for example the quark-antiquark-gluon component of that wave 
function. As we have already pointed out in I we do not expect to see such higher Fock 
states in the decompositions of the Compton amplitude which we have performed. We 
emphasise that they are not missing in our approach. Gluon emissions from the quark 
and antiquark are in our approach contained both in the rescattering term and in the 
dipole-nucleon cross section. It would be interesting to study in detail whether and 
how a consistent separation of such contributions into parts belonging to the photon 
wave function and to the dipole-nucleon scattering amplitude can be done. Naturally 
one would expect that this separation would require to introduce a factorisation scale. 
We recall that in the usual dipole picture the factorisation of the cross section into the 
square of the photon wave function and a dipole-proton cross section does not involve 
such a factorisation scale. 

We consider it an important task for future work to establish a closer relation of 
our approach with perturbation theory. It should be possible to follow our derivation of 
the dipole picture step by step in a perturbative setting. We have for example checked 
our results concerning the different choices for the polarisation vector for longitudinally 
polarised photons in the abelian gluon model of [38] where the cancellations between 
different terms become very transparent. In the comparison with the perturbative 
framework it would be especially interesting to see how the picture of higher Fock 
states of the photon is related to our results. 

Another very important problem is to understand the relation between inclusive 
and diffractive photon-nucleon scattering. This problem has also been addressed in the 
framework of the dipole picture, see for example [31]. In our opinion it is a nontrivial 
question whether the same approximations and assumptions leading to the dipole model 
in inclusive scattering also lead to the usual description of diffractive scattering in 
terms of the dipole-nucleon cross section as given by (147). This question becomes 
particularly challenging when one wants to include corrections from higher Fock states 
in the photon wave function. We expect that the study of diffractive scattering in the 
approach presented here can provide more insight into these problems. 
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We hope that our approach can also help to address other important questions 
concerning high energy lepton-nucleon scattering. We have pointed out for example 
that a study of polarised structure functions in the framework of the dipole picture 
might allow one to investigate and to quantify contributions that are subleading at 
high energies in unpolarised scattering. Another interesting issue is the choice of the 
energy variable in the dipole-proton cross section. Our results strongly suggest that 
the correct energy variable should be the total energy ^/s rather than Bjorken-x. We 
cannot exclude, however, that different approaches to deriving the dipole picture in 
a nonperturbative framework lead to a different result concerning the energy variable. 
We consider this to be an interesting question for future studies. A third problem which 
we consider to be important is the use of the mass-shell condition for quarks. Motivated 
by local parton-hadron duality one often treats quarks produced into the final state as 
on-shell particles, especially in perturbative calculations. The assumption of a quark 
mass shell is also important for the usual formulation of high energy scattering in a 
Hamiltonian approach where the cross sections of higher Fock states of the photon are 
only well-defined for on-shell partons. Due to confinement, however, a mass-shell for 
quarks actually does not exist. We expect that by carefully tracing the effect of using 
the mass-shell-condition for quarks in our framework one can learn something about 
the validity of this assumption. 

The approach presented here might also offer the possibility to study in another way 
the transition region from high to low photon virtualities which has been a key question 
in many applications of the dipole picture to deep inelastic scattering. Most of those 
studies try to use the framework of the dipole picture and introduce modifications of 
the dipole-nucleon cross section at low virtualities. The formulae derived here appear to 
be suited also for an implementation of nonperturbative Green's functions as obtained 
for instance in studies of nonperturbative Dyson-Schwinger equations. In this way one 
could hope to approach the transition region also from the side of low virtualities. In 
practice that would certainly require a number of additional approximations, and would 
hence be less rigorous. Nevertheless, we think that it would be very important from a 
phenomenological point of view to study the transition region between small and large 
photon virtualities using all available methods. 

The general methods developed here can also be applied to other current induced 
reactions at high energies, including exclusive reactions. An interesting application 
would be the study of reactions induced by electroweak currents. The perturbative wave 
functions of W and Z bosons have recently been discussed in [39]. We expect that the 
derivation of the dipole picture for electroweak currents in a nonperturbative framework 
can be performed using essentially the same approximations and assumptions that were 
needed for photon induced reactions. 

Let us finally give a brief summary of the main findings of the present paper and 
of the companion paper I. 

• In I we used functional methods to classify the contributions to the Compton 
amplitude in terms of nonperturbative quark-skeleton diagrams, see figure 1.2. 
The important diagrams at high energies were found to be those of figures 1. 2a 
and I.2b. 

• Cutting the diagrams (a) and (b) of figure 1.2 next to the incoming photon vertex 
and inserting suitable factors of one we could identify parts corresponding to the 
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photon vertex and to quark-antiquark scattering off the proton. Proper account 
was taken of renormahsation and this led to the introduction of the rescattering 
terms, see section 1.2.4. 

• In the high-energy hmit as studied in the present paper the Compton amph- 
tude was spht into a genuine quark-antiquark-proton scattering amphtude and a 
photon-wave-function piece including rescattering terms. Here the assumption of 
the existence of a quark mass shell was made. 

• The equivalence of an incoming or outgoing high-energy photon to incoming or 
outgoing quark-antiquark-dipole states was given explicitly in section 5 for arbi- 
trary reactions. As always, vertex and rescattering terms occurred. 

• The lowest order perturbative wave function of the photon was studied in detail. 
It was shown that the polarisation vectors of the virtual photon have to be chosen 
such that they do not contain terms growing indefinitely as |q| -^ oo. Gauge 
invariance was discussed in detail and was found to be crucial in this context, see 
section 4. 

• In section 6 deep inelastic scattering was treated and the assumptions needed 
to arrive at the usual formulae of the dipole model were spelled out in detail. 
A bound on ai/aT which followed from these assumptions was derived and its 
consequences were discussed in section 7. This bound can be used to test the 
kinematic region in which the usual dipole model can be applied. 

The dipole picture has been surprisingly successful in describing the data obtained 
at HERA and has substantially contributed to the understanding of photon-nucleon 
scattering. In the near future the LHC will offer the possibility to study photon- 
proton and photon-photon collisions at even higher energies for quasi-real photons. 
Particularly suitable for this purpose are highly peripheral heavy ion collisions where 
the Coulomb field of a fast heavy ion acts as a beam of quasi-real photons. We expect 
that the dipole picture will be a very useful tool for analysing the data to come from 
the LHC in this field. A thorough study of the dipole formalism, its nonperturbative 
foundations and its necessary modifications is therefore very important. With the 
present study we hope to have made some useful steps towards this ambitious goal. 
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A Kinematical relations in the high energy Umit 

In this appendix we supplement section 2 with a more detailed discussion of the kine- 
matics relevant to the pinch condition occurring there and further provide some formu- 
lae which are used in the calculation of the photon wave function in section 4. 

A.l The pinch condition 

We start with the discussion of the pinch condition AE — > 0, see (7)- (15). The condition 
(12) on rh?{a,\iT) restricts, for given Q^, the range in a and ky where AE fulfills 
(15). For massless quarks, niq = 0, the condition (12) can be fulfilled for all a with 
< a < 1 with the k^^-range limited as shown in figure 2. For massive quarks, rUg ^ 0, 
the condition (12) can only be fulfilled for a limited range in a. Indeed, for |q| ^ rriq 
this restriction is (as can be seen from (A. 12)) 



OL\ < a < 0.2 



(A.l) 



where 



ai,2 



It 



Arrig 



(A.2) 



This range in a results from requiring that k|, = becomes possible, and the numbers 
given in the following will be calculated for that case. (For a given k^ > to be 
possible the range in a becomes smaller, and its limits are obtained if in (A.2) rni is 
replaced by (k^ + mi).) Clearly, for Q^ = Anii the allowed range in a shrinks to the 
point a = 1/2, whereas for Q"^ ^ co the whole range < a < 1 is allowed. 

In tables 1 and 2 we list some values of ai and 02 for the c- and the 6-quark, 
respectively. We choose as in section 6.5 



rric = 1.3 GeV, 
nib = 4.6 GeV , 



(A.3) 



which implies Q"^ > Am^ = 6.76 GeV^ for the c-quark and Q^ > 4m^ = 84.6 GeV^ for 
the 6-quark. We see from the tables that an allowed a-range (A.l) from 0.05 to 0.95 
requires Q^ > 35.6 GeV for the c-quark and Q^ > 445 GeV for the 6-quark. 



g2[GeV2] 


ai 


02 


20 


0.093 


0.907 


35.6 


0.050 


0.950 


100 


0.017 


0.983 


1000 


0.002 


0.998 



Table 1: Values for Q'^, ai and Q2 for the c-quark 

Finally we want to discuss the actual values of AE in (7) or (10) that one gets in 
various kinematic domains. For high energies we have in the proton rest system 



q° + k° + k'° ^ 2g° = 2v 



Q2 2m^v Q2 



m. 



Q^ 



rripXBj 



(A.4) 
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Q2[GeV2] 


ai 


02 


100 


0.30 


0.70 


445 


0.05 


0.95 


1000 


0.02 


0.98 


10000 


0.002 


0.998 



Table 2: Values for Q^, ai and Q2 for the 6-quark 



where xbj is Bjorken's scaling variable. Inserting this in (15) gives 



< XBj 



rrir 



Q' 



(A.5) 



Let us take now as an example niassless quarks and Q^ = 1 GeV , Q^ = 4 GeV . Then 
AE/rUp < 0.1 for xbj < 0.02. For c-quarks with Q^ = 4GeV2 and 0^ = SGGeV^ we 
get AE/nip < 0.1 for xbj < 0.01. For 6-quarks finally we find with Q^ = 20GeV^ and 
q2 = 450 GeV^ that AE/rUp < 0.1 is only reached for xbj < 0.005. 

With these numbers for xbj we do by no means want to establish limits for the 
validity of the dipole picture. We only want to show how small AE is in concrete cases. 



A. 2 Further kinematical relations 

In the high energy limit |q| -^ oo and with the conditions (36) we find from (4)-(6) the 
expansions 

I k'i + mf, ,_ / (k|, + m'^y ^ 



A;° = a|q| + - 



Vj, 



2 a|q| 



+ 



and 



k'" = (l-a)\q\ + 



1 kl + mi 



2 (l-a)|q| 
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m 



(K+m^y 



We further have the products 



fc°fc'° = a(l - a)|qp + i ( -^ + '-^) (k^ + m^) + O ( iii±^ 



2 Vl-a 



a 



(A.6) 



(A.7) 



{A.i 



and 

(/c° + mg)(/c'° + mq) = a(l-a)|qp+mg|q| + - ( :;-^ + ^^) (k| + m2)+m^ 
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(A.9) 



The leading contribution to {AE) "^ (see (3) and (7)) at high energy is obtained as 

gO + fcO + fe'O 



(A^)-^ 



a(l-a)(gO + |q|) 
a(l - a)Q'^ + k|, + m2 



l + O 



k^ + m^ 



(A.IO) 
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and we have 
Furthermore, we find 



2|q| 

_ k'^ + m'^ 
a(l — a) 



l + O 



k% + ml 



Vj, 



(A.ll) 



(A.12) 



for which < a < 1 is required, as was assumed throughout the discussion of the high 
energy hmit, see (9) and (36). 

For the Fourier transformation of the photon wave function from transverse mo- 
mentum space to coordinate space we need the well-known integral 



d'^h 



(2vr) 



T JkyRy 



^^ = ^KoMt). 



(A.13) 



The computation of the wave function of transversely polarised photons also requires 
the integral 

which is obtained from (A.13) by differentiating with respect to the components of Ry. 

B Separate gauge invariance of different contributions to 
the Compton amplitude 

In this appendix we discuss the separate gauge invariance of different contributions to 
the skeleton decomposition of the Compton amplitude, see figure 1.2. Here we deal in 
particular with those contributions that are subleading at high energies. 

We first recall from I that the different contributions in the decomposition (1.13) of 
the Compton amplitude can be expressed as (see (I.A.9)-(I.A.17)) 
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where j = a, . . . ,g, and the explicit expressions for the J'^^' are given in appendix LA. 
For p and p' off shell we can integrate in (B.l) by parts with respect to y and y' and 
get 
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Here, as usual, one has to take the hmit fl -^ nip and ^ -^ nip for the expression in 
square brackets first and then to multiply with the Dirac spinors. The contraction of 
J\^U)t^'^ with Qiy can then be obtained in analogy to (20) as 



quM^^r{p',P,q) 



2TTnipZp 
- ^ + nip) 



Us' (pO (- / + rup) AA(^') (/, p, q) 



nip 
m„ 



Usip) , 



(B.3) 



where we have defined 



M^^\p',p,q) 
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jii) 
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(B.4) 



From (B.3) we see that a nonvanishing contribution to q^MS^'^^ is obtained only if 
A^") has poles both at ^' = nip and at |( = nip. 

In section 3.1 we have already discussed the separate gauge invariance of A^''*' and 
M^'. Let us now consider the gauge invariance of the amplitude M.^^'. Following 
steps similar to those in (21) we find that 
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As a consequence we obtain from (LA. 16) with (LA. 11) 
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and inserting this into (B.4) we arrive at 
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We can now interpret the integrals over y' and y as Fourier transformations. The first 
term in the curly brackets in (B.7) leads to a contribution corresponding to an incoming 
proton of momentum p + q and an outgoing proton of momentum p'. Consequently, 
this contribution has poles at ^' = nip and at ^+ /^ = nip, but not at ^ = nip. 
Similarly, the second term in the curly brackets gives rise to a contribution which has 
poles at ^ = nip and ^'— ^ = nip, but not at ^' = nip. Therefore, neither of the two 
contributions can give a nonvanishing contribution when inserted in (B.3), and hence 



quMiPf''{p',p,q)=0. 
This completes the proof that A^ '-•'-* is separately gauge invariant. 



(B. 
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In order to discuss the remaining parts of the amphtude it wih be useful to consider 
first the u- and the d-quark contributions to the electromagnetic current, 



J^{x) = Quu{xh\{x), 
4{x) = Qdd{x)-i^d{x) , 



(B.9) 
(B.IO) 



and their matrix elements between incoming and outgoing proton states. Using the 
LSZ reduction formula and the methods explained in I we get in analogy to (I- 10) for 
these matrix elements 
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In analogy to the treatment of (I. A. 8) in appendix LA we obtain using (LA. 7) and 
(LA. 11) for the functional integrals in (B.ll) involving these currents 
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We now turn to the amplitudes M^'^' and Ai^'^^ which are related to J'^'^' and J'^'^\ 
respectively, via (B.l). From (LA. 13) we find using (B.5) 
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(B.14) 
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where £ is given by 
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Again using (B.5) we find from (LA. 14) tliat also the derivative of J^'^' with respect to 
x^ can be expressed in terms of iS in a similar way, 
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We see that in the sum of (B.14) and (B.16) the terms with the delta function involving 
x' cancel. Hence we have 



._d_ 

' dx^ 



j{c) ^ j{d) 



5(4)(a._y)_<j{4)(y/_^y 



£. 



(B.17) 



From (B.4) we then find using (B.12) and (B.13) 



Qd( (My')Jdio)i^piy)) +Qd( V'p(y')v^p(y)Tr (-4^)(o,o;G)7 

' G,q,q \ \« / / G 



+ Q.((V'p(y')^^(o)V'p(y))^^_ + Q«(^p(y')^p(y)Trfi4")(o,0;G)7^ 



(B.18) 



With this we have expressed the sum of M^'^' and AA*-"^' in terms of the proton field 
-i/jp. Hence the four terms in square brackets can be understood as contributions with 
incoming and outgoing protons, so that we can apply the same argument as in the case 
oiAf^-'' above. We again interpret the integrals over y' and y as Fourier transformations 
and find that none of the two contributions arising from the two terms in curly brackets 
has simultaneous poles at ^' = nip and j5 = nip. Consequently, M^'^' + M^'^' gives a 
vanishing contribution when inserted in (B.3), such that the sum of M^'^' and Ai^ ' is 
gauge invariant by itself. 



qjM^;l^''{p',p,q)+M'^friP'^P^l) 



0. 



(B.19) 



Note that neither M^'^' nor M^'^' are separately gauge invariant. If we consider for 
example Ai^^' we see that after inserting (B.14) into (B.4) the delta function containing 
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x' will give rise to a term that has simultaneous poles at -^ = nip and ^ = nip and hence 
will contribute to q^M^'^'^'^. Only due to the cancellation of this contribution with the 
opposite one in q^M^'^'^^'^ is the sum of A^*-"^^ and Ai^'^' separately gauge invariant. 

Next we consider the gauge invariance of the amplitude TW*-^-*. From (LA. 15) we 
find with (B.5) 

. d 



' dx'^ 



J 



(e) 



6^^\x - y) - 6^*\y' - x) 



^ ex' (3' ^''^ af3^ 



(d) 



Qu{Qu + Qd)-S';>Jy',y;G) 



^^sPiy',x';Gh^-/;\x',y;G)^ ^ ^S^^^^^iy' ,y;G) 
- (a' ^ /?') - (a ^ /?) + {a' ^p',a^ (5) 



-^Fc'a^y'-' y^ G)-Sy^,p{y', y; G) - {a ^ (3) 



(B.20) 



G 






G,q,q 



+ QJ Vp(y')V'p(y)Tr ( -SP{0,0;G)r 



Inserting this in (B.4) and expressing it in terms of the proton field ipp via (B.12) and 
(B.13) we obtain 



(B.21) 



G 

V ^ ' <Ji'?i'? \ V * J I G 

Again in complete analogy to the discussion of }yv-'' above we find that this expression 
does not contain terms which have simultaneous poles at fl = nip and p = mp, and 
hence does not contribute to q^M^^'^'^ when inserted in (B.3), 

g.A^2'''(p',P,9) = 0, (B.22) 

so that A^'*^-* is separately gauge invariant. 

Finally we turn to the amplitude MS^'. Here inspection of (LA. 17) shows that the 
x-dependence of J^'-^' and hence of MS^> is fully contained in the same factor 



r 



j;Q,(-l)Tr(7-i4«)(x,x;G)) 



(B.23) 



which occurred already in J^^> . We can therefore in analogy to (24)-(26) derive that 

g,A^if^"(/,p,(?)=0. (B.24) 

Thus also MS^> is separately gauge invariant. 
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